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ABSTRACT 


The cost of operating ships is difficult to predict. A historic ship's 
operating cost database is maintained by the Miltary Sealift Command 
(MSC); but, it is very difficult to extract or manipulate the data to support 
prediction or regression analysis. An alternative was sought that would 
reduce the effort for the user when attempting to make predictions from the 
data. If the data for each cost category (salary, training, fuel, port and 
miscellaneous, subsistence, ship's equipage, and voyage repairs) could be 
well approximated using probability distributions, then the costs of an 
operational scenario, with estimates of the uncertainties, could be obtained 
through use of a Monte Carlo simulation. 

The MSC data was divided into two subsets, one for model fitting and 
one for validation. Once probability distributions had been fit to the data, a 
Monte Carlo simulation tool was developed using the Crystal Ball® 
simulation add in to Microsoft Excel®. The data analysis and cost model 
were then validated using the empirical data. 

Based on the results, the Cost Simulation model provides a useful tool 
for predicting operating costs and supports sensitivity analysis of various 


ship's operating cost scenarios. 
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THESIS DISCLAIMER 

The reader is cautioned that computer programs developed in this 
research may not have been exercised for all cases of interest. While every 
effort has been made, within the time available, to ensure that the programs 
are free of computational and logic errors, they cannot be considered 


validated. Any application of these programs without additional verification 


is at the risk of the user. 
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EXECUTIVE SUMMARY 


The cost of operating ships is difficult to predict. A historic ship's 
operating cost database is maintained by the Military Sealift Command 
(MSC); but, it is very difficult to extract or manipulate the data to support 
prediction or regression analysis. An alternative was sought that would 
reduce the effort for the user when attempting to make predictions from the 
data. If the data for each cost category (salary, training, fuel, port and 
miscellaneous, subsistence, ship's equipage, and voyage repairs) could be 
well approximated using probability .distributions, then the costs of an 
operational scenario, with estimates of the uncertainties, could be obtained 
through use of a Monte Carlo simulation. 

The MSC data was divided into two subsets, one for model fitting and 
one for validation. The model fitting subset was analyzed using graphical 
data analysis. The object of this effort was to fit probability distributions to 
the data. Good probability distribution fits were found for each of the cost 
categories examined in the subset. 

Once probability distributions had been fit to the data, a Monte Carlo 
simulation tool was developed using the Crystal Ball® simulation add-in to 
Microsoft Excel®. The simulation tool was designed with a user interface to 


reduce the technical knowledge needed by the user to operate the application. 


The data analysis and cost model were then vahdated using the 
empirical data. The simulation tool was run and the results compared to the 
actual per diem rates in use at MSC. This analysis showed that the 
simulation model produced results close to the actual per diem rate. The 
results pointed to either possible inaccuracies in the overhead rates used at 
MSC, or problems with the simulation model. 

The simulation model was run without indirect (overhead) costs 
included and compared to the historical direct costs in the entire database. 
These runs showed the simulation model to accurately predict direct costs. 

On the basis of the results, the Cost Simulation model provides a 
useful tool for predicting direct operating costs and supports sensitivity 
analysis of various ship's operating cost scenarios. Further study is required 
in the area of indirect (overhead) costs to permit use of the simulation model 
for prediction of total costs of operation. This would enable MSC to use the 


simulation tool for setting per diem rates with a higher degree of accuracy. 





I. BACKGROUND AND PROBLEM STATEMENT 
Military Sealift Command, United States Pacific Fleet (MSCPAC) is a 
civilian manned and operated shipping line owned by the United States 
Government. The mission of MSCPAC is fourfold: 

1. Operate and maintain Naval Fleet Auxiliary Force (NFAF) ships to 
provide direct support to U.S. Navy combatant ships. 

2. Operate and maintain Prepositioned Sealift ships as required by the 
National Command Authority (NCA). 

3. Operate and maintain Special Missions Force ships as required for 
specialized military purposes such as oceanographic and hydrographic 
surveys, undersea surveillance, and missile telemetry collection and range 
instrumentation. 

4. Charter and contract ocean cargo services as necessary to support the 
military commitments. 

The first three of the missions are performed primarily by civilian 
crewed ships that are owned and operated by the government. To accomplish 
this, MSCPAC operates a fleet of thirty-four ships at an annual cost of $500 
million. 

In recent years, the funding scheme was changed from the Industrial 
Fund to the Defense Business Operating Fund (DBOF). Under DBOF 
funding, the sponsor for each ship is billed for the services provided by the 


] 


ship. The charges are billed on a cost per day or "per diem" basis. The 
sponsor pays these bills using his Operating and Maintenance, Navy Funds 
(OMN). MSCPAC must meet all of 1ts operating costs using funds obtained 
in this manner. 

In the post cold-war atmosphere of dwindling defense budgets, the 
sponsors are in some cases demanding a lower per-diem charge. In the case 
of Special Mission ships, the sponsor of the cable ships, Navy Space 
Surveillance and Warfare Command (SPAWAR), has threatened to contract 
with a civilian contractor who 1s bidding a lower price. 

Any overcharges reduce the precious OMN funds that are needed to 
conduct all the sponsor's missions. Eventually, when the ship sponsor's 
budget is cut further, he will be forced to fulfill his requirements with the 
lowest bid, which frequently is the civilian contractor. This increased need 
for cost competitiveness requires the ability to forecast budget needs to a 
greater degree of accuracy than is presently common in government practice. 
It is clearly in the best interest of MSCPAC to preserve their viability as an 
economical alternative to the civilian contractors by obtaining the ability to 
make tighter budget forecasts. The costs of operation be analyzed and cut to 
the minimum to permit MSCPAC to quote the lowest per-diem rate possible. 

There are three major categories of cost in the MSCPAC operation, the 
cost of overhead for facilities and shore based employees, the cost of 
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chartering and contracting cargo services, and the cost of operating the 
government owned and operated ships. In the opinion of the MSCPAC 
comptroller and staff, the major difficulties are 1n the third category. 

The first category should be fairly predictable and should not contribute 
greatly to the total cost of MSCPAC. 

The second category comprises a different form of control by MSCPAC. 
During the negotiating of a charter or contract for services, the best price 1s 
sought. When finally billed for the charters and contracts, MSCPAC pays 
the predetermined price. MSCPAC is later reimbursed by the sponsor for the 
services provided by MSCPAC through the charter or contract. This area 
normally does not present a problem to the accounting department. 

The final category, the cost of operating the government owned and 
operated ships is the area of most concern to the comptroller. 

But with existing data and tools, questions from management for cost 
analysis information become painful forays into voluminous heaps of data 
that are difficult to interpret and analyze. Often, the answers require 
several man days effort to collate and assemble the data into a manageable 
form. The results are not timely enough to meet the demands for information 
on which to base decisions. 

Presently, all payment vouchers written by MSC are paid by MSC 
Headquarters in Washington, DC. The transactions are recorded on a 
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mainframe based accounting system called the Financial Management 
Information System (FMIS). The area headquarters such as MSCPAC can 
access data from FMIS through a personal computer landline linkup system 
called FMIS Gateway. 

The PC/mainframe interface is cumbersome. The information is 
accessible on-line, but it is very difficult to analyze one screen at a time. All 
reports are routinely downloaded from FMIS Gateway in the form of 
printouts. When the need for analysis arises, the accounting department 
must gather vast amounts of paper records. Frequently, by the time an 
answer is produced, the question is no longer of interest! 

Management needs a more responsive tool for budgeting and conducting 
"what if" analysis. This tool should aid in budget forecasting by allowing 
management to predict the effects of anticipated fiscal changes with a 
reasonable amount of accuracy. The ability to conduct sensitivity analysis 
would permit the tightening of costs in areas that would most affect the total 
cost of operation. 

Another advantage of a tool of this type would be the ability to analyze 
year to date expenditures. If certain fiscal policies can be predicted to 
outspend the budget before the fact, those practices can be changed to allow 
the activity to remain within budget. At present, this ability is almost 
non-existent other than the traditional "stubby pencil" methods. 
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The purpose of this thesis is to develop and demonstrate a cost 
simulation tool through which this process can be streamlined so answers 
can be found within a matter of minutes rather than days. The tool is based 
on the premise that each of the categories of operating cost for the 
government owned and operated ships behaves in a way that can be modeled 
using known probability distributions. Once these probability distributions 
are known, cost studies for different periods of times and conditions can be 
made very quickly using a personal computer and Monte Carlo simulation. 

The probability and simulation models are combined in a popular 
spreadsheet program with an add-in simulation solver. This powerful 
combination brings the ability to make decisions based on solid forecasting 
directly to the manager's desktop. No longer will the manager be totally at 
the mercy of the accounting department when he needs rapid answers to his 
questions. 


The remainder of this thesis 1s arranged as follows: 


1. Chapter II explains the methodology used in the project. 

2. Chapter III gives a description of the components of each cost category. 

3. Chapter IV presents the data analysis. The chapter is rather lengthy 
and is summarized in the first section. A more detailed reading of the data 
analysis is the subject of the remainder of the chapter. 

4 Chapter V explains the cost analysis model. 

5.Chapter VI briefly explains the Monte Carlo simulation method as 


applied in this cost simulation model. 


6. Chapter VII presents a validation of the data analysis. 
7. Chapter VIII gives the conclusions. 
8. The appendix gives the visual basic code used in building the enhanced 


user interface for the spreadsheet model. 


П. METHODOLOGY 

Exploratory data analysis and Monte Carlo simulation are combined to 
build a cost simulation tool to predict the direct costs of operating the MSC 
owned and operated ships. Monte Carlo simulation is a scheme in which 
random numbers are used for solving certain stochastic or deterministic 
problems where the passage of time plays no substantive role. In this case, 
random numbers from a uniform distribution on the unit interval, U(0,1), 
will be processed through probability distribution transformations to model 
the operating cost category data according to relationships found by 
performing data analysis on the historical operating cost data. 

The Monte Carlo simulation model is implemented as an add-in to the 
spreadsheet commonly used at MSCPAC, Microsoft Excel®. The Monte Carlo 
simulation is implemented in Crystal Ball® from Decisioneering Inc., 
Denver, Colorado, which runs as an add-in to Microsoft Excel. This system 
runs on any PC that is capable of running Microsoft Excel. In most cases at 
MSCPAC, the software will be run on a 386 PC. 

This approach requires an extensive data analysis of historical cost 
data. The data is in the form of monthly expense records and is analyzed to 


obtain probability distributions to model each category for each ship class. 


The central assumptions behind the cost simulation tool are that the 
ships operating costs are random in nature and can be modeled by known 
probability distributions. 

A spreadsheet database is first created in Microsoft Excel to enable the 
transfer of data to the data analysis program, A Graphical Statistical System 
(AGSS). A Graphical Statistical System (AGSS) is used at the Naval 
Postgraduate School under a test site agreement with IBM Research. We are 
indebted to Dr. Peter Welch for making this possible. 

Due to software incompatibility between the FMIS Gateway system and 
the PC spreadsheets, the data was hand keyed into the database. Future 
enhancement of the cost simulation tool could be achieved by establishing 
direct links between FMIS Gateway and the spreadsheet database. 

The actual data analysis is accomplished in AGSS. Each ship class cost 
category becomes a variable in AGSS. "These variables are first analyzed 
using descriptive plots such as histograms and kernel density estimates to 
determine which theoretical distributions might be appropriate. The 
variables are next plotted and analyzed using the probability distribution 
fitting capability of AGSS. Each variable is plotted in various manners and 
goodness-of-fit statistics are calculated to determine the best probability 
distribution fit. Figure 1 depicts the three plot view used for probability 


distribution fitting. 
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Figure 1 
Sample Probability Distribution Plot 

In the upper left hand plot, the data 1s plotted as a histogram of relative 
frequency of occurrence versus cost. This histogram is superimposed upon a 
plot of the selected probability density which has been generated by AGSS. 
AGSS computes the maximum likelihood estimates of the distribution 
parameters, together with asymptotic confidence regions. 

The upper right plot shows the data set plotted as an empirical 
cumulative frequency distribution superimposed on the cumulative 


distribution function for the fitted probability distribution. In this view, 


AGSS offers the ability to plot Kolmogorov-Smirnov bounds as a further aid 
in determining a good fit between the data and the theoretical distribution. 

The lower left hand plot in Figure 1 shows a probability plot. In this 
plot, the theoretical distribution 1s represented as a straight line diagonally 
across the plot. The data 1s plotted and superimposed such that if 1t was a 
perfect fit to the theoretical distribution, the plot would fall exactly on the 
line. 

The lower nght hand section of Figure 1 contains tabular data that 
helps determine goodness of fit. Four different goodness of fit tests are 
provided for uncensored data. These are the Chi-Square, 
Kolmogorov-Smirnov, Cramer-von Mises, and Anderson-Darling tests. Due 
to the size of the figures, the tabular data for probability distribution plots 
used in this analysis is found in Appendix A. 

By analyzng the plots and comparing them with those of other 
distributional fits, itis possible to determine the best distributional fit for the 
data. These distributional fits are used as distributional assumptions in the 
simulation model. 

The cost simulation model is written in the form of a Microsoft Excel 
workbook. Each worksheet page in the workbook is used for a different ship 
class or model variation. The ship classes examined in this project will be 
the T-AO 187 class fleet oiler and T-ATF 166 class oceangoing tug. 
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Microsoft Excel is used primarily for the convenience of the intended 
users. MSCPAC personnel are familiar with Microsoft Excel so use of Excel 
for this model will reduce the training requirements for the users. 

Microsoft Excel supports the use of the Crystal Ball Monte Carlo 
simulation add-in. This product is available in versions for Microsoft Excel 
and Lotus 1,2,3. Crystal Ball is a user-friendly, graphically-oriented 
forecasting and risk analysis program. Through Monte Carlo simulation, 
Crystal Ball forecasts the entire range of results possible for a given 
situation. Crystal Ball uses probability distributions to describe the 
uncertainty in the assumption cells of the spreadsheet model. There are 
sixteen different probability distributions available in Crystal Ball for 
describing the relationships 1n the model. 

Microsoft Excel allows the user to operate in the Microsoft Windows 
environment. The Windows graphical user interface (GUI) allows users to 
quickly grasp complex concepts in a more intuitive manner. All planned cost 
simulation tool users are already skilled 1n using Microsoft Windows based 
programs. 

Microsoft Excel also supports a rich macro language (Visual Basic) to 
allow automating the application. This macro language is used to develop a 
simplified interface for the users. The primary goals for the cost simulation 
tool are accurate results and ease of use. 
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Pairwise correlation between each cost category is examined using the 
Spearman Rank correlation coefficient. Where there 1s any noted pairwise 
correlation, it is compensated for in the assumption parameters of Crystal 
Ball. Multiple correlation is also examined. 

The final facet in establishing the accuracy and creditability of the cost 
sumulation tool is to thoroughly validate the model and the distributional 
assumptions. The first method of validation is to divide the data and analyze 
each subset independently to determine whether distributional assumptions 
are valid for all subsets of the data set. In this case, the first half of the data 
set (chronologically divided) is used to obtain initial assumptions and the 
second chronological half of the data is used to validate the assumptions. 

Finally, test run simulations are conducted to determine the accuracy of 
the cost simulation that corresponds to the actual conditions for randomly 
selected ships and months. The results of the simulations are compared to 


the actual month end account figures. 
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III. MILITARY SEALIFT COMMAND COSTS 

MSC accounts for ship operating costs using direct costs and indirect 
costs. The direct costs are made up of cost categories that are economically 
feasible to trace to a particular ship (or cost center). Indirect costs are 
overhead costs that are not economically feasible or practical to trace to an 
MSC ship. For example, the cost of xerox copy paper for use in the 
headquarters offices are paid for by DBOF revenues from the ships but its 
consumption cannot be directly traced to a particular ship. 

Indirect costs are difficult to accurately predict. Indirect costs from 
previous years are examined and a budget estimate is made based on the 
amounts spent in previous years. An indirect cost budget for the coming 
fiscal year is then allocated as a percentage of estimated direct costs for the 


year. The total cost charged to a sponsor is simply: 


Total Cost = Direct Costs + Allocated Indirect Costs 


This analysis focused only on the direct costs in an effort to show a 
relationship between the historical database and known probability 
distributions. This permits the use of Monte Carlo simulation for predicting 


direct costs. 


13 


The overhead costs (designated as administrative indirect costs for this 
analysis) were not precisely known by MSC and were not directly related to 
the actual operation of the ships. The Cost Simulation Tool cannot simulate 
these costs since one of its assumptions 1s that costs are directly related to 
the operation of the ships. The overhead rate used in the model was given by 
MSC. No attempt was made by this analysis to simulate these costs. 

Overhead costs are allocated as a percentage of direct cost throughout 
the Department of Defense’. This method is also used at MSC. An overhead 
rate of 19% of the direct cost was.given by the MSCPAC Director of 
Operations as the current rate being used by MSC. The calculation of the 


rate is as follows: 


Overhead Rate = (Budgeted Indirect Costs+Estimated Direct Costs)x 100 


A. INDIRECT COSTS. 
Three categories of indirect cost are allocated to two overhead cost 
pools. The first pool of overhead costs used wil be known as the 


administrative pool and includes the following: 


i Headquarters and other administrative support costs. 
2. Physical plant and building costs. 


‘Young, Douglas, "Complexities, Impact of Overhead", U.S. Army Comptroller 
Office, Resource Management, Ist Quarter, 1994 
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This pool makes up about 65% of the model's budgeted indirect cost or 
about 19% of the actual direct costs. This pool consists all of what would be 
referred to as overhead costs by MSC. 

The second pool of indirect costs used in the model consists of planned 
maintenance and docking costs. This pool is used in the Cost Simulation 
Tool model but in actual practice these costs are budgeted direct costs 
because they are traceable to a particular ship and are budgeted in advance. 
The model assumption that a ship is 100% available to the sponsor conflicts 
with the actual case (i.e., that the ships are not available to the sponsor while 
undergoing planned maintenance and/or drydocking.) The probability 
distribution modelling feature of the Cost Simulation Tool is only accurate 
for cases in which the ship is 100% available to the sponsor, therefore the 
maintenance and docking costs will be treated as if they are an indirect cost. 
The planned maintenance costs are budgeted in advance and included in the 
estimate of total direct cost for the coming year. The estimated budgeted 
overhead rate (administrative indirect cost pool in the model) and the 
estimate of total direct cost (which includes the planned maintenance and 
drydocking costs) are used in determining the per diem rate that is charged 
to the sponsors for use of the ships. Since MSC lumps the costs together in 


their estimate for per diem rate determination, it is reasonable for this 
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analysis to model the maintenance and drydocking costs as budgeted indirect 
costs which will be summed to the direct costs in the model. 

No effort was made to obtain probability distributions to fit the 
planned maintenance and drydocking costs. The historical maintenance 
costs were summed and allocated as an average percent of total direct cost 
for each class of ships for the entire period of time for which the data was 
obtained. For each model run, this allocated percentage is added to the 
percentage used for the other cost pool. The maintenance and docking costs 
make up about 35% of the model's budgeted indirect costs and 10.4% of the 
total direct costs. Since these costs are planned for and known in advance, 
they are not treated as a random variable. This slight departure from actual 
practice will permit us to better use the model to predict the overall operating 
costs. 

The budgeted indirect costs will be allocated using percentage of total 
direct cost as the basis. In actual MSC accounting practice, the overhead 


costs are also allocated as a percentage of total direct cost. 


B. SALARY COSTS. 
This category includes all payroll costs incurred in paying seagoing 
employees. The following is a detailed list of salary costs as listed on actual 


FMIS printouts: 
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. Salary data communication. 


Civilian mariner base pay. 

Civilian mariner overtime pay. 

Civilian mariner premium pay. 

Civilian mariner hazardous duty pay. 

Civilian mariner beneficial suggestion awards. 
Civilian mariner incentive pay. 

Civilian mariner awaiting assignment pay. 
Civilian mariner indoctrination pay. 


. Military pay. 

. Civilian mariner annual leave earned. 
. Civilian mariner sick pay earned. 

. Civilian mariner shore leave earned. 

. Civilian mariner health insurance. 

. Civilian mariner life insurance. 

. Civilian mariner retirement fund. 

. Civiian mariner FICA. 


C. TRAINING COSTS. 


This category includes damage control and safety training costs 


obtained in route to the ship by officers and crew. The following is a detailed 


list of training costs: 
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. Officers damage control instructor school. 
. Civilian mariner firefighting training enroute. 


Civilian mariner small arms training. 
Civilian mariner safety training. 


. Maritime academy cadet training. 
. Miscellaneous training. 
. Closed circuit television system training. 


D. FUEL COSTS. 


This category includes the cost of diesel fuel and petroleum lubricants 


and greases. 
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E. SUBSISTENCE COSTS. 
This category includes the cost of subsistence for civilian officers and 
crew. Military officer and crew subsistence is subsidized by the military pay 


account and is not a cost to MSC. 


F. PORT AND MISCELLANEOUS COSTS. 
This category includes the cost a number of items. The following is a 


detailed listing of port and miscellaneous costs: 


1. Transportation of items to and from ship. 
2. Consumables. i 

3. Spare parts. 

4. ADP supplies. 

5. Software. 

6. Medical supplies. 

7. ADP equipment. 

8. Docking and other fees. 

9. Piloting and towage. 

10. Panama canal tolls. 

11. Utilities. 

12. Security guards. 

13. Civilan mariner repatriation travel. 
14. Civilian mariner other travel. 

15. Laundry. 

16. Movies/tapes. 

17. INMARSAT. 

18. Medical expenses. 

19. Other miscellaneous expenses. 
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G. SHIP'S EQUIPAGE 
This category includes the cost of equipage items such as binoculars, 


tools, and foul weather clothing. 


H. VOYAGE REPAIRS 

This category includes the cost of unplanned repairs completed during 
the course of normal ship operation. This item was extracted from the 
category of maintenance and repair costs (which consists of planned 
maintenance, drydocking and voyage repairs). The remainder of the 
maintenance and repair items ( planned maintenance and drydocking) are 


included in the budgeted indirect cost category. 
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IV. DATA ANALYSIS 

A significant portion of the work behind the cost simulation tool 
lies in data analysis. Financial cost data is maintained on a mainframe 
computer located at Military Sealift Command (MSC) headquarters in 
Washington, DC. This detailed information includes the disposition of every 
payment voucher written by MSC. The information is accessible in read only 
form by MSCPAC and the other area commands via a Personal Computer 
(PC) modem hook up using software called the Financial Management 
Information System (FMIS) PC Gateway. 

Data for this project was obtained from FMIS for all MSCPAC 
ships operated during fiscal years 1992 and 1993. The data consists of 
monthly summaries for the various cost categories required to operate ships. 
The total cost of operating a MSCPAC ship can be regarded as the sum of 


the following categories: 


Salaries. 

Training. 

Fuel and lubricants. 

. Subsistence. 

Port and miscellaneous (including spare parts). 
Ship's equipage. 

Maintenance and repair. 

. Budgeted overhead. 
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The monthly totals for each cost category were entered into a 
Microsoft Excel O spreadsheet database. This spreadsheet database allowed 
for ease in manipulating data to form variables to be used in the graphical 
data analysis software, A Graphical Statistical System © (AGSS). 
Additionally, Microsoft Excel is presently in use at MSCPAC, so every effort 
was made to stick to its use to reduce the amount of training needed for the 
future users of this application. 

To ease the process of filtering the data to remove various 
undesirable data elements, Microsoft Excel's database query feature was 
used. For example, if it was desired to remove all data from months in which 
the ship was not available to the sponsor for 100 percent of the time, the 
query would be made requesting that those cases be deleted. The raw data 
can be similarly manipulated for any other desired case. 

Graphical data analysis allows the characteristics of the data to 
be studied visually as well as through computational statistics. Sometimes it 
is far easier for the eye to see a relationship than it is to discern it from 
computational results. The cost category data can be modeled with known 
probability distributions and this relationship is easy to see graphically. 

AGSS! is an interactive system for both two and three 
dimensional scientific-engineering graphics, applied statistics, and data 


‘IBM Corporation; A Graphical Statistical System (AGSS) : An 
Introduction; pg 1. 
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analysis. | AGSS allows the user to create graphics, explore data 
interactively, analyze data using functions of applied statistics, develop 
customized graphics functions, and manage, review, and store work sessions. 
In this case, AGSS was chosen for its excellent capability in the area of 
distribution fitting. 

AGSS has capabilities for fitting any of 18 univariate probability 
distributions to a set of data. The system computes maximum likelihood 
estimates as well as several other estimates of the distribution parameters, 
together with asymptotic confidence regions, and it can produce three 
dimensional contour and surface plots of the likelihood function. 

AGSS provides a number of graphical comparisons of the 
empirical data with the theoretical fitted distribution. These displays help 
the user judge visually how well different distributional assumptions apply 
to the data. They provide graphs on which some representation of the 
theoretical distribution is superimposed on the corresponding empirical plot. 
Examples of this are plots of the histogram and fitted density function, 
empirical and fitted cumulative distribution functions (CDF), and probability 
plots. In a probability plot, if the empirical data corresponds exactly to the 
quantiles of the theoretical distribution, the points wil lie exactly on the line 


that runs diagonally across the plot from lower left to upper right. 
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Kolmogorov-Smirnov (KS) bounds, at any desired confidence level, can be 
superimposed on the CDF and probability plots. 

Four different goodness of fit tests are provided for uncensored 
data. These are the Chi-Square, Kolmogorov-Smirnov (KS), Cramer-Von 
Mises (C-VM), and Anderson-Darling (AD) tests. These tests give a 
quantitative measure of the goodness of fit. These values can be compared to 
the visual interpretation of the fit to help verify the quality of the fit. In 
graphical data analysis though, the visual impressions often give the analyst 
a far better indication of a probability distribution fit than do these statistical 
tests. 

Figure 2 is a probability distribution plot created by AGSS. This 
view illustrates many of the principles discussed in the following paragraphs. 
Three plots are shown in this view. The plot in the upper left hand corner is 
a superimposed plot with a histogram of the empirical data superimposed on 
a plot of the theoretical probability density function selected for the fit. The 
plot in the upper right corner shows a plot of the empirical CDF 
superimposed on the theoretical CDF of the distribution selected for the fit. 
The plot in the lower left hand corner is a probability plot where the 
empirical percentiles are plotted against corresponding percentiles of the 
theoretical probability distribution. In this plot, if the empirical data exactly 
matched the theoretical probability distribution quantiles, all of the points 


23 


Т-АО 187 Class Salary 
LOGISTIC DENSITY FUNCTION, N=25 LOGISTIC CUMULATIVE DISTRIBUTION FUNCTION, N=25 


е 
ч“ 


0.4 0.8 0.8 


CUMULATIVE PROBABILITY 


0.2 





0 
0 


0.4 0.5 1.2 0.4 0.5 1.2 
Dollars x 964.076 Dollars x 964,075 


LOGISTIC PROBABILITY PLOT, N=25 





0.4 0.5 12 
Dollars x 964,075 


Figure 2 
Sample Probability Distribution Plot 


would plot on the straight line. A table of statistical and goodness of fit 
results normally is shown in the lower nght hand corner of the plot. This 
four way plot is used extensively in the data analysis. 

KS bounds are included on the CDF and probability plots. A 
theoretical CDF or straight line on the probability plot passing outside the 
bounds indicate lack of fit. The KS values are calculated for the user 
specified confidence level (95 percent in this analysis). The KS test statistic 
is based on the maximum difference between the observed empirical CDF 


and a hypothesized distribution across all values of x. As KS statistic values 
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decrease, the fit between the sample and the theoretical distribution is 
improved. 

The significance level output by AGSS is the p-value used in 
statistics. The significance level is the smallest level of significance of the 
test statistic for which the null hypothesis would be rejected. In the cases 
examined in this project, the null hypothesis is that the data came from the 
distribution specified for the fit. The higher the significance level the more 
likely that the fit will be accepted. 

The Cramér-Von Mises (C-VM) test statistic is based on the 
integral of the squared distance between the empirical and theoretical 
curves. The C-VM value, like the KS value, should be as low as possible. 
The C-VM significance level is given only in ranges for the C-VM statistic. 

The Anderson-Darling (AD) test statistic is an attempt to 
overcome a drawback in both the KS and C-VM tests. Both the KS and the 
C-VM are not sensitive to departures from the null hypotheses that occur in 
the tails of the distribution. This is improved in the Anderson-Darling test 
by using a weighted distance measure, the weight being the reciprocal of the 
standard deviation of the difference between the curve functions. The 


smaller the AD value, the more likely the fit is accepted. 


? Lewis, Orav; Simulation Methodology for Statisticians, Operations 


Analysts, and Engineers: Volume 1; pg 369 
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Note that the Chi-Square test results are reported by AGSS for 
some of the probability distribution fit plots. The Chi-Square test is used 
when fitting discrete distributions or when continuous distributions can be 
modeled by grouping the data into mutually exclusive discrete bins such as 
in a histogram. The normally recognized lower limit for frequency 1n a bin 
for the Chi Square test is five. With the small size of the data sets (25 and 16 
data points), this limit is tested 1n nearly every case, and the Chi-Square test 
will therefore not be used in this analysis. 

The remainder of this chapter 1s organized as follows: 

A. Assumptions - Description of the assumptions used in 

analyzing the MSCPAC data. 

B. Summary of Analysis Results. 

C. T-AO 187 Class Data Analysis - a detailed account of the data 

analysis of each cost category for ships of the T-AO 187 class. 

D. T-ATF 166 Class Data Analysis - a detailed account of the 


data analysis of each cost category for ships of the T-ATF 166 class. 


A. ASSUMPTIONS. 

For each data set (or cost category for a particular class of ships) 
numerous plots were made in an attempt to find which of several probability 
distributions best fit the data. Early in this process, it became clear that 
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most of the cost categories were subject to a fair degree of variability, 
depending on the operating schedule for the ship. 

For example, a ship that 1s deployed performing services for the 
sponsor burns significantly more fuel than a ship of the same class that is in 
port frequently due to maintenance and administrative requirements. The 
same sort of variability holds in the other six categories of direct cost 
considered in the model (salaries, training, subsistence, port and 
miscellaneous costs, ship's equipage, and voyage repairs). 

Fortunately, the MSCPAC Operations Department recently 
conducted a study? in which the operating tempos of MSCPAC ships were 
examined. In this study, it was decided to consider a ship in one of three 
states: available to sponsor, not available to sponsor due to maintenance 
requirements, and not available to sponsor due to administrative or training 
requirements. The state for each of the assigned ships was recorded for the 
last two fiscal years, and projected for the next five fiscal years. 

When assigning a per-diem rate for the ships, one must assume 
that the ship is exclusively available to the sponsor for the period which is 
being paid for (i.e., ships available to sponsor 100 percent of the time he is 
paying for) After making this assumption, the data sets were scrutinized 


and only the months in which the ship 1s 100 percent available to the sponsor 


° Military Sealift Command; MSCPAC Operations Department QMB, 
October 1993 
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(1.е., no days spent in satisfying maintenance or administrative 
requirements) were included in the data set. This filtering resulted in a 
significant reduction in variance of the observed costs. 

Since the purpose of the model is to predict cost, 1t is reasonable 
to assume that in an operational status, no cost category should show a 
negative or zero total for a month. This assumption is important because the 
data set contains some instances in which zero and negative sums are carried 
as a matter of convenience for the accountants who later shift funds between 
accounts. This assumption permits filtering the data set to remove any 
months in which a negative or zero balance is shown for any category. The 
filtering removes some points that may bias the plots toward the Nu 
data. 

Rather than including the entire maintenance and repair 
category in the direct cost portion of the model, only the maintenance and 
repair account for voyage repairs is included as a direct cost. This is done 
since the other accounts under maintenance and repair are not normally 
used in months where the ship is 100 percent available to the sponsor. Those 
accounts are maintenance and dry-docking costs which, for our purposes, are 
summed over the entire period and then divided to permit allocation of the 


costs for the period of time concerned in the model run. The voyage repair 
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account, on the contrary, takes into account the normal voyage repairs that 
are charged to the ship while it is in operation. 

For this project, only the T-AO 187 class oilers and the T-ATF 166 
class oceangoing tugs are analyzed. In both cases, there are a sufficient 
number of ships of these classes assigned to MSCPAC to comprise data sets 
of a large enough size that data analysis and the inferences drawn from the 
data analysis will apply. 

Each ship class 1s analyzed separately. For these ship classes, the 


following cost categories are analyzed: 


. Salaries. 

Training. 

. Fuel and lubricants. 

. Subsistence. 

Port and miscellaneous. 
Ship's equipage. 

. Voyage repairs. 
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B. SUMMARY OF RESULTS 
The theoretical distributions that were selected as the best 
approximations for each category are given in Table 1 for the T-AO 187 class 
and Table 2 for the T-ATF 166 class. The maximum likelihood estimates for 
the parameters, mean, standard deviation, and a subjective assessment of 


the quality of the fit are also given. In some cases it was difficult to choose 
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from a small set of distributions, and in other cases, none of the distributions 
seemed to fit very well. In every case, the distribution chosen is the one the 
author judged to best explain the data when all evidence was considered. 
Both data sets were analyzed to determine if there was any 
correlation between cost categories. In only one pair of categories, T-ATF 166 
Training versus Subsistence costs, is a significant correlation found. This 


correlation is incorporated into the Cost Simulation Tool model. 


Table 1. 
T-AO 187 CLASS DATA ANALYSIS RESULTS 


Cost Category Distribution Parameters Mean Std Dev Fit Quality 
Salary Logistic о:=692,601, В =95,974 692,601 174,179 Good 
Training Weibull С =1.5615, а=24,864 22,346 14,622 Good 
Fuel Logistic 0.—202,026, B —53,822 202,026 97,624 Fair 
Subsistence Logistic о=5,966, В =3,723 19,486 4,944 Marginal 
Port and Misc. Weibull С =2.7352, а=171,010 152,570 61,177 Excellent 


Ship's Equipage Weibull С =1.0616, а=21,037 20,414 18,890 Good 
Voyage Repairs Gamma 0:—0.838, B 2127,700 106,950 97,542 Good 


Table 2. 
T-ATF 166 CLASS DATA ANALYSIS RESULTS 


Cost Category Distribution Parameters Mean Std Dev Fit Quality 


Salary Normal и=123,880, © =54,864 123,880 54,864 Good 
Training Gamma 0.21.3843, B 222243 3,104 2,689 Marginal 
Fuel Weibull С=1.7563, а =49,847 44,409 26,105 | Good 
Subsistence Logistic =3,433,8=795.7 3,670 1,894 Marginal 
Port and Misc. Lognormal =10.075,0=0.6869 30,047 23,333 Marginal 
Ship's Equipage Lognormal р-=7.4952,6=13561 4,513 10,879 Good 


Voyage Repairs  Lognormal ц=11.082, с =1.068 115,010 167,790 Good 
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For a detailed coverage of the analysis of each cost category for 
each class of ship, the reader 1s encouraged to read the following sections of 
this chapter. Following the cost category analyses, the correlation analysis 1s 
found for each class of ship. 

C. T-AO 187 CLASS ANALYSIS 

1. Salary Cost. 

Figure 3 shows the probability distribution fit plot for 


T-AO 187 class salary cost. The plot shows a good fit using the logistic 
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Figure 3 
T-AO 187 Class Salary Cost Probability Distribution Plot 
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distribution. The theoretical CDF and distribution fall well within the 
chosen Kolmogorov-Smirnov (KS) bounds. 

The goodness of fit tests indicate a good fit with KS value 
of 0.11668 and KS significance level of 0.88546. For the salary data, the 
C-VM value of 0.63255 and significance level > 0.15 indicate a good fit. For 
this plot, the AD value of 0.4681 and significance level of > 0.15 indicate that 
this is an adequate fit. 

Although the fit looks good, there is some departure in 
the tails. This is seen in the plots of Figure 3. The values sharply start in 
the left tail and also sharply drop off in the right tail as seen in both the CDF 
view and the probability plot. Figure 4 is a view in which the empirical 
density has been superimposed over the theoretical density (empirical is solid 
curve, theoretical is dash-dot curve). In the superimposed plot, the departure 
in the tails is also seen as well as some difference in the area before the peak 
of the curve. 

Figure 5 is a multiple box plot of the salary data. The 
actual data box closely resembles that of the small random sample box plot. 
The small random sample box plot was plotted from a random sample (using 
the fitted distribution and parameters) of the same size as the actual data 
set. The large random sample box plot was made from a random data set ten 
times the size of the actual data set. The large random sample box plot 
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Figure 4 
T-AO 187 Class Salary Cost Superimposed Density Plot 
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Figure 5 
T-AO 187 Class Salary Cost Multiple Box Plot 
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closely resembles the actual data set plot except that the tails are longer due 
to the increased size of the data set. 

The long lines extending above and below the actual data 
set box indicate some symmetry in the tails, although the lines are of slightly 
different length. The horizontal line across the box indicates the median and 
in this case since it 1s so far above the mean it shows some lack of symmetry 
near the peak of the curve. This 1s also seen in Figure 4. The data 1s skewed 
slightly to the left as indicated by the mean circle lying below the median 
line. 

Figure 6 1s a symmetry plot that agrees with the above 
analysis. The points le below the у-х line which indicates that the data is 
skewed to the left. 

The fit with the logistic distribution is the best obtained 
for the salary cost data. The normal distribution is the closest runner up, but 
the goodness of fit results are not as strong. The fitted parameters for the 
logistic distribution are $692,601 and $95,974 for alpha and beta, 
respectively. The mean and standard deviation to be used in the Cost 
Simulation Tool Model are $692,601 and $174,179 respectively. 

2. Training Cost. 

Figure 7 is the probability distribution fit plot for monthly 

training cost data. The histogram seems to be almost tailor-made for the 
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Figure 6 
T-AO 187 Class Salary Cost Symmetry Plot 
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Weibull density function. The default bin selection by AGSS resulted in a 
histogram that 1s very similar to the shape of the density. Both the CDF and 
probability plots show the theoretical density to lie well within the KS 
bounds. The points in the probability plot of Figure 7 lie very close to the 
hne. 

The KS value of 0.09619 with significance of 0.97484 and 
the C-VM value of 0.031256 with significance > 0.15 indicate an excellent fit 
for the data with the Weibull distribution. The AD value of 0.24434 with 
significance > 0.15 indicates that the fit is good. 

Figure 8 is the superimposed density plot showing the 
empirical data as a solid curve and the theoretical curve as a dot-dashed 
curve. Some lack of fit in the tails is seen in this view. The second peak seen 
in the empirical curve is due to one high leverage point in the data at 
approximately $45,000. 

Figure 9 shows a multiple box plot for the data. The nght 
skewness of the data set can be clearly seen in this view. This would be 
expected for data that were sampled from the Weibull distribution. The left 
tail is steeper and shorter than the right tail as indicated by the shorter line 
below the box. The median line is below the mean and shifted toward the 
bottom of the box. The small random sample box is slightly higher than the 
actual data set and large random sample data set boxes. Since the actual 
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Figure 8 
T-AO 187 Class Training Cost Superimposed Plot 
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Figure 9 
T-AO 187 Class Training Cost Multiple Box Plot 
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data set resembles the large random data set, the fit is good. 

Figure 10 is a symmetry plot that also indicates a right 
skewness of the data. This is seen by the fact that the distances are much 
higher to points above the mean than to points below the mean. 

The Weibull distribution will be used as the distribution 
for T-AO 187 class monthly training cost in the Cost Simulation Tool model. 
The gamma distribution is a close runner up, but, the Weibull fits are much 


better. The shape parameter to be used in the model will be 1.5615 and the 
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Figure 10 
T-AO 187 Class Training Cost Symmetry Plot 
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scale parameter will be $24,864. The fitted mean and standard deviation for 
the data were $22,346 and $14,622, respectively. 
3. Fuel Cost. 

Figure 11 1s a probability distribution fit plot for the fuel 
cost data. The best fit for the data is obtained with the logistic distribution. 
Although the superimposed density histogram plot shows what appears to be 
a good fit, it must be remembered that the appearance of the histogram can 
be significantly altered by changing the number of bins used. In this case, 
AGSS default values for bin selection provide a histogram that is a good 
match for the density. The CDF and probability plots indicate a lack of fit in 
the tails. The theoretical distribution is fully contained within the KS 
bounds for the data. The KS value is 0.15724 with a significance level of 
0.56679. The C-VM value is 0.10111 with significance level > 0.15. The AD 
value of 0.71454 with significance level > 0.15 indicates an adequate fit. 

Figure 12 is the superimposed density plot for the fuel 
cost data. As indicated above, the peak and middle area of the empirical 
density (solid curve) fit the theoretical density (dot-dashed curve) very well, 
but, the tails of the empirical density are thicker than those of the theoretical 
density. 

Figure 13 is the multiple box plot for the fuel cost data. It 


shows that the majority of the empirical data set is tighter than either of the 
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Figure 11 
T-AO 187 Class Fuel Cost Probability Distribution Plot 


random sample data sets. The view shows less similarity between the 
empirical data and the random samples than is indicated in Figures 11 or 12. 

The spread of the actual data box is similar in both 
directions with outhers on both the left and right ends of the curve. Figure 
14, the symmetry plot for the fuel cost data, also indicates strong symmetry 
with the points lying relatively even about the y=x line. 

The logistic distribution will be used bi the Cost 
Simulation Tool model to represent fuel cost data. There were no other 


distributions that exhibited a good fit for the data. The fitted parameters are 
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Figure 12 
T-AO 187 Class Fuel Cost Superimposed Density Plot 


$202,026 and $53,822 for alpha and beta, respectively. The fitted mean and 


standard deviation values are $202,026 and $97,624, respectively. 


4. Subsistence Cost. 
Figure 15 is the Probability Distribution Fit plot for 


subsistence cost. The superimposed histogram and density plot shows that 
the majority of the data lies 1n the two center bins. The CDF and probability 
plots show that the fit 1s less than perfect, but the theoretical distribution 


does Пе within the KS bounds. The trace from the data departs from the 


theoretical shapes 1n the center due to the concentration of data there. The 
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Figure 13 
T-AO 187 Class Fuel Cost Multiple Box Plot 
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Figure 14 
T-AO 187 Class Fuel Cost Symmetry Plot 
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Figure 15 
T-AO 187 Subsistence Cost Probability Distribution Plot 


KS value of 0.21565 and significance level of 0.19535 indicates that there is a 
fit for the logistic distribution, albeit not a very strong one. Both the C-VM 
and AD suggest that there is an adequate fit, but the values for both tests, 
0.27196 and 1.5123 respectively, are relatively large when compared with 
other values in this analysis. In both cases, the significance levels are > 0.15. 

Figure 16, the superimposed density plot, indicates a 
much better fit than was evident in the previous plot, but the tails of the 


empirical density are much thicker than those of the theoretical density. 
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Figure 16 
T-AO 187 Subsistence Cost Superimposed Density Plot 


Even more of the probability is concentrated in the center of the empirical 
plot than is characteristic of the logistic distribution. 

Figure 17 is the multiple box plot for the subsistence 
data. The actual data box does not remotely resemble the boxes for the 
random sample cases. The box and the extended lines of the actual data box 
are compressed and clearly indicate that the majority of the data is in the 
narrow peak of the distribution. There are five points recognized as distant 


outliers (the solid circles seen above and below the box). 
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Figure 17 
T-AO 187 Class SubsistenceCost Multiple Box Plot 


Figure 18 is a symmetry plot for the data set. The plot 
indicates a symmetry which is consistent with the fit for the logistic 
distribution. This was also seen in the superimposed density plots. 

The logistic distribution fit for this data is not very well 
supported, but it is the only fit obtainable from the 18 distributions available 
in AGSS. The logistic distribution will be used in the Cost Simulation Tool 


Model. The fitted parameters are $5,966 and $3,723 for alpha and beta. The 
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Figure 18 
T-AO 187 Class Subsistence Cost Symmetry Plot 


fitted mean and standard deviation were $19,486 and $4,944, respectively, 


and will be used as parameters in the model. 


9. Port and Miscellaneous Cost. 

Figure 19 is the Probability Distribution Fit plot for port 
and miscellaneous cost. The three plots all show a strong relationship 
between the data and the Weibull distribution. The CDF and probability 
plots show that the theoretical distribution is wholly included within the KS 


bounds. The stepped data traces very nearly mimic the CDF. The data 
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Figure 19 
T-AO 187 Port and Misc.Cost Probability Distribution Plot 


points all plot very close to the y—x line of the probability plot. The KS value 
of 0.089326 with significance level of 0.98844 indicate an excellent fit for the 
Weibull distribution. The C-VM value of 0.026164 and significance level of 
> 0.15 also indicate a strong fit. The AD value of 0.1888 with significance 
level of >0.15 also indicate an adequate fit. 

Figure 20 is the superimposed density plot for the data. 
The empirical plot 1s very close to that of the theoretical density. A slight 


lack of fit in the tails 1s visible 1n this view. 
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Figure 20 
T-AO 187 Port and Misc.Cost Superimposed Density Plot 


Figure 21, the multiple box plot for the port and 
miscellaneous cost, indicates a fairly symmetric distribution in the actual 
data. The actual data plot is similar to the large random sample plot, 
indicating a fit exists. The small random sample plot is slightly less 
symmetric and has a higher mean than the actual data plot. Figure 22 is the 
symmetry plot for the data. The points lie mostly below the y=x line, but they 
are fairly close to the line which accounts for the symmetry. 

The fit of the Weibull distribution to the data set is a 
strong one and will be used in the Cost Simulation Tool Model. The fitted 


shape and scale parameters are 2.7352 and $171,010, respectively. The 
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Figure 21 
T-AO 187 Port and Misc. Cost Multiple Box Plot 
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Figure 22 
T-AO 187 Port and Misc. Cost Symmetry Plot 
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mean and standard deviation are $152,570 and $61,177, respectively. 


6. Ship's Equipage Costs. 

Figure 23 1s the Probability Distribution Fit plot for ship's 
equipage cost. The data exhibits a good fit for the Weibull distribution. The 
CDF and probability plots clearly show that the theoretical distribution lies 
within the KS bounds. The stepped empirical CDF tracks fairly close to the 


theoretical CDF. The data points are close to the line in the probability plot. 
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Figure 23 
T-AO 187 Ship's Equipage Cost Probability Distribution Plot 
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The KS value of 0.10514 with a KS significance level of 0.94511 indicates a 
good fit exists. The C-VM value of 0.044199 with significance level > 0.15 
also indicate a good fit. The AD results indicate an adequate fit exists with 
AD value of 0.29121 and significance level of > 0.15. 

The superimposed density plot shown in Figure 24 
displays less of a fit between the empirical density and the theoretical 
density than the test statistics would lead one to believe. There is clearly a 


lack of fit in both the left and right tails. 
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T-AO 187 Ship's Equipage Cost Superimposed Density Plot 
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Figure 25 1s a multiple box plot for the ship's equipage 
cost data. The data appears to be skewed to the right as indicated by the 
shorter line below the box and the fact that the median is near the bottom of 
the box. The actual data plot closely resembles the large random sample 
plot. The small random sample plot 1s slightly different than the other two. 

Figure 26 is a symmetry plot for the data set. It also 
supports the right skewness of the data with the long right tail indicated by 
the points lying above the line. The distance to points above the mean is 
much greater than the distance to points below the mean. 
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Figure 25 
T-AO 187 Ship's Equipage Cost Multiple Box Plot 
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Figure 26 
T-AO 187 Ship's Equipage Cost Symmetry Plot 


The fit for the Weibull distribution will be used for this 
category in the Cost Simulation Tool Model. There is also a good fit for the 
gamma distribution with this data, but the fit is not as strong. The fitted 
shape and scale parameters of 1.0616 and 21,037, respectively, for the 
Weibull distribution will be used in the model. The mean and standard 


deviation for the data were $20,414 and $18,890, respectively. 
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7. Voyage Repair Cost. 

Figure 27 is the probability distribution fit plot for the 
voyage repair cost data. The superimposed histogram and density plot is not 
of much use in this case due to lack of detail in the default view presented by 
AGSS. The relative frequency scale 1n this view seems to be too large for the 
data presented by the histograms. A good fit with the gamma distribution 1s 
indicated by the information in this figure. The CDF and probability plots 
both show that the theoretical distribution is contained within the KS 
bounds. The KS value of 0.11317 with significance level of 0.90803 is 
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Figure 27 
T-AO 187 Voyage Repair Cost Probability Distribution Plot 
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another indicator of a good fit. The C-VM value of 0.072187 and significance 
level > 0.15 also indicates a fair fit. The AD indicates an adequate fit with a 
value of 0.47551 and significance level > 0.15 respectively. 

The superimposed density plot of Figure 28 shows less of 
a resemblance between the empirical density and the theoretical density 
than the test results above seem to indicate. The lack of fit in the tails is also 
clearly seen in this view. The data also is skewed to the nght in this view. 

Figure 29 is a multiple box plot of the voyage repair cost 


data. The right skewness of the actual data is very pronounced in this view. 
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Figure 28 
T-AO 187 Voyage Repair Cost Superimposed Density Plot 
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Figure 29 
T-AO 187 Voyage Repair Cost Multiple Box Plot 


The left tail is short as evidenced by the short line extending from the bottom 
of the box. The right tail, on the contrary, is very long as seen by the long 
line extending from the top of the box. The median line is toward the bottom 
of the box and lies below the mean. The random sample views are not of 
much help in this case. Based on this plot, an inference that a fit does not 
exist would be drawn. 

Figure 30 is a symmetry plot that also shows the long 
right tail by the difference in the axes and the fact that the points plot above 


the y=x line and depart further from the line as the cost increases. 
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Figure 30 
T-AO 187 Voyage Repair Cost Symmetry Plot 


The gamma distribution fit will be used for voyage repair 
cost in the Cost Simulation Tool Model. The runner up is a fit with the 
lognormal distribution that is not as strong as the fit with the gamma 
distribution. The fitted parameters for the gamma distribution are 0.83755 
for alpha and 127,700 for beta. The mean and standard deviation are 
$106,950 and $97,542, respectively. 

8. Correlation Between Data Categories. 
The multiple correlation test feature in AGSS is used to 


test for correlation between any pair of the categories of data. The presence 
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of any correlation would be important to explain during the cost analysis 
phase of the project. If correlations exist, they will also require 
compensation in the Cost Simulation Tool Model. Crystal Ball, the 
simulation add-in software for Microsoft Excel, allows the user to input 
Spearman Rank correlation between any two data sets used 1n a model. 

The multiple correlation test for the T-AO 187 data 
matrix indicates that no significant correlation exists between any pair of 


data sets in the data matrix. Table 3 is the multiple correlation matrix for 


the T-AO 187 Class data. 


Table 3. 
T-AO 187 CLASS DATA MULTIPLE CORRELATION MATRIX 


Salaries Training Fuel Subsistence Port Equipage Voyage Repairs 
Salaries 1 0.0694 0.191 0.246 | 0.339 -0.251 0.137 
Training 0.0694 1 -0.177 -0.121 0.055 -0.228 0.138 
Fuel 0.191  -0.177 1 -0. 118 0.142 -0.172 -0. 13 
Subsistence 0.246  -0.121  -0.118 1 -0121 -0056 0.096 
Port 0.339 0.055 0.142 -0. 121 1 -0149 0.146 
Equipage -0.251 -0.228 -0.172 -0.056 -0.149 1 0.153 
Voyage Repairs 0.137 0.138 -0.13 0.096 0.146 0.158 1 


D. T-ATF 166 CLASS ANALYSIS. 


1. Salary Cost. 


Figure 31 is the probability distribution fit plot for the 


T-ATF 166 class salary cost data. A normal distribution fit was obtained for 


58 


T-ATF 166 Class Salary Cost 
NORMAL DENSITY FUNCTION, N=16 NORMAL CUMULATIVE DISTRIBUTION FUNCTION, N=16 


e 
= 


0.4 0.6 
0.8 


0.8 


0.2 
0.4 


RELATIVE FREQUENCY 
0.3 
CUMULATIVE PROBABILITY 


0.2 


о о 


Dollars 


NORMAL PROBABILITY PLOT, N-16 








Figure 31 
T-ATF 166 Class Salary Cost Probability Distribution Plot 


this data set. The superimposed histogram-density plot shows a wide 
dispersion of values with an exceptionally large density in the center bin of 
the histogram. The CDF and probability plots show that the theoretical 
distribution is completely contained within the KS bounds. The KS value of 
0.15001 and significance level of 0.86422 indicate a fair fit between the 
empirical data and the theoretical density. The C-VM results of 0.075901 
and significance > 0.15 also support this hypothesis. The AD test indicates a 


fit exists with a value of 0.47923 and significance > 0.15. 
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Figure 32 1s the superimposed density plot for the salary 
cost data. As mentioned above, the plot shows that the fit with the normal 


distribution is generally good, but that there is a lack of fit in the tails of the 
distribution. 


Figure 33 is the multiple box plot for the data. This view 
indicates that the actual data is similar to the large random sample, but 


shghtly different in spread from the small random sample box. 


Figure 34 is the symmetry plot for the salary cost data. A 


right skewness is indicated by the information in this plot. The distances to 
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Figure 32 
T-ATF 166 Class Salary Cost Superimposed Density Plot 
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Figure 34 
T-ATF 166 Class Salary Cost Symmetry Plot 
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points below the median are generally greater than the distances to points 
above the median. 

The normal distribution fit will be used in the Cost 
Simulation Tool model. There are no other distributions that are considered 
good fits for this data set. The fitted mean and standard deviation are 
$123,880 and $54,864, respectively. 

2. Training Cost. 

Figure 35 is the probability distribution fit plot for the 
training cost data. An acceptable fit was obtained using the gamma 
distribution. The superimposed histogram-density plot shows that the 
histogram has the same general shape as the gamma density. The CDF and 
probability plots show that the theoretical distribution is contained within 
the KS bounds. The stepped empirical CDF does not adhere closely to the 
theoretical CDF in some areas of the plot. The same is true in the probability 
plot where the data points cross back and forth over the line in the right tail 
of the plot. The KS value of 0.20696 with significance level of 0.4996 does 
not give as strong a case for a fit as some of the other data sets seen earlier. 
The C-VM value of 0.08997 with significance level >0.15 supports the 
hypothesis that a fit exists. The AD value of 0.4847 with significance level of 


> 0.15 indicates that a fit exists. 
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Figure 35 
T-ATF 166 Class Training Cost Probability Distribution Plot 


Figure 36 1s the superimposed density plot for the data. 
In this view, the left part of the data set up to the peak seems to be the area 
where there 1s a departure from the theoretical gamma density. 

Figure 37 1s a multiple box plot which shows the data to 
be right skewed. This agrees with what was seen in the superimposed 
density plot. The data would be left skewed if not for the two largest values 
in the data set. Thisis shown by both of the random sample views. 

Figure 38 1s a symmetry plot that also shows the strong 


right skewness of the data. The distance to most of the points above the 
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T-ATF 166 Class Training Cost Superimposed Density Plot 
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Figure 37 
T-ATF 166 Class Training Cost Multiple Box Plot 
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Figure 38 
T-ATF 166 Class Training Cost Symmetry Plot 


median is less than $1000, whereas the distances to points below the median 
ranges from $1250 to $2500. 

The gamma distribution fit will be used in the Cost 
Simulation Tool model for T-ATF 166 monthly training cost. A fit is also 
obtained using the Weibull distribution, but the fit is not as good as this fit 
with the gamma distnbution. The fitted parameters for the gamma 
distribution are 1.3843 for alpha and 2,242.50 for beta. The mean and 


standard deviation were $3,104.40 and $2,638.50, respectively. 
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3. Fuel Cost. 

Figure 39 is the probability distribution fit plot for the 
fuel cost data. These views indicate a good fit of the fuel cost data to the 
Weibull distribution. The superimposed histogram-density plot shows that 
there is a good resemblance between the histogram and the theoretical 
Weibull distribution. The stepped empirical CDF tracks fairly closely to the 
theoretical CDF. The data points plot close to the line 1n the probability plot. 
The theoretical distribution is contained within the KS bounds in both the 
CDF and probability plots. The KS value of 0.13903 and significance of 
0.91656 support the hypothesis that a fit exists. The C-VM value of 0.032114 
and significance level of > 0.15 also support the hypothesis. The AD value of 
0.211 with significance level of > 0.15 indicate a fit exists. 

Figure 40 is the superimposed density plot. The shape of 
the empirical density is very close to that of the theoretical density with some 
small departures near the peak and in the right tail. Both the empirical and 
theoretical densities are skewed slightly to the nght. 

Figure 41 is the multiple box plot of the data. There is 
similarity between the actual data box and the random sample boxes. This 
view agrees with the slight right skewness seen in the superimposed density 


plot earlier. The line extending from the bottom of the box is smaller than 
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T-ATF 166 Class Fuel Cost Probability Distribution Plot 
the line extending above the box indicating a longer right tail. The median is 
near the middle of the box and is very close to the mean value. 

Figure 42 is the symmetry plot of the data. The points 
below the median are all closer to it than the points above the median. This 
also agrees with the slight right skewness and long right tail shape discussed 
above. 

The Weibull distribution fit will be used for fuel cost data 
in the Cost Simulation Tool Model. An acceptable fit is also obtained using 


the gamma distribution, but the fit is not as good as the Weibull distribution 
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Figure 40 
T-ATF 166 Class Fuel Cost Superimposed Density Plot 


fit. The fitted parameters for the Weibull distribution were 1.7563 for shape 
and 49,847 for scale. The mean and standard deviation were $44,409 and 
$26,105, respectively. 

4. Subsistence Cost. 

Figure 43 is the probability distribution fit plot for the 
subsistence cost data. The plots indicate that a fair fit exists for the data 
with the logistic distribution. The superimposed histogram-density plot 
shows that the histogram shape 1s similar to that of the density. The stepped 


empirical CDF departs somewhat from the theoretical CDF due to the high 
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T-ATF 166 Class Fuel Cost Multiple Box Plot 
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Figure 42 
T-ATF 166 Class Fuel Cost Symmetry Plot 
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Figure 43 
T-ATF 166 Class Subsistence Cost Probability Distribution Plot 


density of data in the center, but the theoretical CDF is completely contained 
within the KS bounds. In the probability plot, the points are not very close to 
the line although the line is completely contained within the area of the KS 
bounds. The KS value of 0.23775 with significance of 0.32626 does support 
the hypothesis that a fit exists. The C-VM value of 0.2226 and significance 
level > 0.15 do not rule out a fit. The AD value of 1.2804 with significance 
level > 0.15 indicate that a fit exists. 

Figure 44 is the superimposed density plot for the 


subsistence cost data. This view presents a much stronger case for using the 
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Figure 44 
T-ATF 166 Class Subsistence Cost Superimposed Density Plot 
logistic distribution fit. 


The empirical density very nearly duplicates the 


shape of the theoretical density. There is some slight departure in the tails 
of the distribution. 


Figure 45 is the multiple box plot for the data. This plot 


shows that symmetry exists for the data set. The actual data is more tightly 
squeezed in the middle of the distribution than the random sample boxes. 
The tight grouping of the data does conform to the tight peaked appearance 


of the logistic distribution though it 1s slightly tighter which accounts for the 


crossing back and forth over the line in the probability plot. 
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Figure 45 
T-ATF 166 Class Subsistence Cost Multiple Box Plot 


Figure 46 is the symmetry plot for the data. This plot 
also indicates symmetry for the data with only one data point departing from 
the vicinity of the y—x line. 

The logistic distribution will be used for the T-ATF 166 
subsistence cost in the Cost Simulation Tool model. There are no other 
distributions that exhibited a good fit for this data. The fitted parameters 
are 3432.7 for alpha and 795.66 for beta. The mean and standard deviation 


are $3669.80 and $1894, respectively. 
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Figure 46 
T-ATF 166 Class Subsistence Cost Symmetry Plot 


5. Port and Miscellaneous Cost. 

Figure 47 is the probability distribution fit plot for port 
and miscellaneous cost. From the analysis in this plot there is a marginal fit 
to the lognormal distribution. Most of the data is grouped tightly about 
$20,000, but there is one value of over $400,000 and two values in the 
$100,000 range that exert leverage to shift the density toward the extreme 
right. A satisfactory fit with any distribution was unattainable until the 
point at $400,000 was removed. A marginal lognormal fit was then obtained. 


With the aforementioned caveat, the theoretical distribution in the CDF and 
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Figure 47 
T-ATF 166 Port and Misc. Cost Probability Distribution Plot 


probability plots falls entirely within the KS bounds. The KS value of 
0.239932 with significance of 0.35671 satisfy the hypothesis that a fit exists. 
The C-VM value of 0.18118 with significance level » 0.15 does not rule out 
the hypothesis that a fit exists. The AD value of 1.0456 with significance 
level > 0.15 also indicate a fit. 

The superimposed density plot shown in Figure 48 indicates a 
lack of fit where there is a gap in the empirical data between the clump of 


values at about $20,000 and the large values mentioned previously. 
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Figure 48 
T-ATF 166 Port and Misc. Cost Superimposed Density Plot 


Figure 49 is a multiple box plot of the port and miscellaneous cost 
data. The distribution is mostly a tight bunch of values as can be seen by the 
tightly squeezed box at the lower values. There are also 4 outlier points 
which account for the long nght tail of the data. The actual data box is 
similar to the large random sample box. 

Figure 50 is a symmetry plot that shows that the points located in 
the tight bunch are symmetric, but points in the right tail cause a radical 


departure from the symmetry. 
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Figure 49 
T-ATF 166 Port and Misc. Cost Multiple Box Plot 


The only fit obtainable, as mentioned before, is achieved after 
removing the point at $400,000. Discussion with both the MSCPAC 
accounting and operations departments indicates that $400,000 for one 
month's port and miscellaneous cost 1s much greater than normal and would 
not be expected to be repeated. The lognormal distribution will be used in 
the Cost Simulation Tool model. The fitted parameter values are 10.075 and 
0.68694 for mu and sigma, respectively. The values for the mean and 


standard deviation are $30,047 and $23,333, respectively. 
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Figure 50 
T-ATF 166 Port and Misc. Cost Symmetry Plot 


6. Ship's Equipage Cost. 

Figure 51 is the probability distribution fit plot for the 
ship's equipage data. This data was found to satisfy a fit for the lognormal 
distribution. The superimposed histogram-density plot shows a good 
resemblance between the histogram and the theoretical density. The stepped 
empirical CDF tracks fairly well with the theoretical CDF. The data points 
fall close to the line in the probability plot. The theoretical distributions lie 
completely within the KS bounds in both plots. The KS value of 0.13769 


with a significance level of 0.92204 indicates that a fair fit exists. The C-VM 
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Figure 51 


T-ATF 166 Ship's Equipage Cost Probability Distribution Plot 


value of 0.050035 with significance level » 0.15 also supports the hypothesis 


that a fit exists. The AD result of 0.31448 with significance > 0.15 indicates a 


Figure 52 is the superimposed density plot which shows that 


there 1s some similanty in the shape of the empirical density and the 


theoretical density. There are a few data points 1n the $10,000 range that 


cause a lack of fit 1n the right tail. 


Figure 53 1s the multiple box plot for the ship's equipage data. 


This plot shows that the data are right skewed with a shorter line below the 
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Figure 52 
T-ATF 166 Ship's Equipage Cost Superimposed Density Plot 


box and the median below the middle of the box as well as being below the 
mean value. The right tail is much longer than the left tail primarily due to 
the data points in the $10,000 range. The actual data box is similar to the 
large random sample box plot. 

Figure 54 is the symmetry plot of the data. This view also 
supports the idea that the data are right skewed. Most of the distances to 
points above the mean are longer than distances to points below the mean. 

The lognormal distribution will be used to model ship's equipage 


cost in the Cost Simulation Tool model. The Weibull distribution also fit the 
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Figure 53 
T-ATF 166 Ship's Equipage Cost Multiple Box Plot 


data, but the fit with the lognormal distribution is much better. The fitted 
parameters for the lognormal distribution are 7.4952 and 1.3561 for mu and 
sigma, respectively. The fitted mean and standard deviation are $4512.80 
and $10,379 respectively. The large standard deviation is cause for concern 
here. The addition of more data points after more operation time will help in 


re-plotting, fitting, and finding a fit with reduced variance. 
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Figure 54 
T-ATF 166 Ship's Equipage Cost Symmetry Plot 


7. Voyage Repair Cost. 

Figure 55 is the probability distribution fit plot for the 
voyage repair cost data. A fit was obtained with the lognormal distribution. 
The superimposed histogram-density plot shows that although there is some 
resemblance between the histogram and the data, there are also gaps in the 
data at several places. The stepped empirical CDF does exhibit a fit with the 
theoretical CDF. The data points plot near the line in the probability plot. 
The theoretical distribution in both plots lies within the KS bounds. The KS 


value of 0.12223 with significance of 0.9706 indicates that a good fit exists. 
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Figure 55 
T-ATF 166 Voyage Repair Cost Probability Distribution Plot 


The C-VM value of 0.052749 with significance level > 0.15 supports the fit. 
The AD value of 0.3337 with significance level > 0.15 indicates a fit exists. 

Figure 56 is the superimposed density plot for the voyage 
repair cost data. The empirical density seems to have a good match on the 
left side of the curve, but on the right side the data points with large gaps in 
between cause a departure from the natural shape of the theoretical 
distribution. 

Figure 57 1s a multiple box plot of the data. The right 


skewed nature of the data is clearly seen in this view where the line below 
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Figure 56 
T-ATF 166 Voyage Repair Cost Superimposed Density Plot 


the box is significantly shorter than the one above the box. This indicates 
that the night tail is much longer than the left tail. The median line is well 
below the middle of the box and lies below the mean as well. The actual data 
box resembles the small random sample except that the actual data box is 
slightly narrower. 

Figure 58 is a symmetry plot of the data. This view 
clearly shows the right skewness of the data. The distance to points above 
the mean is much greater than the distance to points below the mean. All of 


the points plot well above the y=x line. 
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Figure 57 
T-ATF 166 Voyage Repair Cost Multiple Box Plot 
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Figure 58 
T-ATF 166 Voyage Repair Cost Symmetry Plot 
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The lognormal fit for this data set will be used in the Cost 
Simulation Tool model. The data also fit the Weibull distribution, but the fit 
was not as good. The fitted parameters for the lognormal distribution are 
11.082 and 1.068 for mu and sigma. The mean and standard deviation are 
$115,010 and $167,790, respectively. 

8. Multiple Correlation Test Results. 

The T-ATF 166 class data matrix was examined using the 
multiple correlation function in AGSS to determine if there were any 
significant correlations between data set pairs. Table 4 is the pairwise 
correlation matrix. 

From the matrix in Table 4, there are three cases of 
possibly significant correlation. Correlation seems to exist between salary 
cost and voyage repair cost (0.396), training cost and subsistence cost (0.372), 
and port and miscellaneous cost and voyage repair cost (0.639). Each of 
these correlations is further examined in a pairwise manner. In each case, a 
bivariate scatter plot is plotted. Then linear regression is performed to 
permit further graphical examination of the relationship between the data 
sets. In all three cases, there is one high leverage point which seems to be 
driving the correlation. In each case, the high leverage point is removed, the 


data is replotted, and the correlation is recalculated. 
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Table 4 
T-ATF 166 CLASS MULTIPLE CORRELATION DATA MATRIX 


Salary Training Fuel Subsistence Ports & Misc Equipage Voyage Repairs 


Salary 1 0.091  -0.05 -0. 152 0.094 0.0388 0.396 
Training 0.091 1 0.108 0.372 -0.204 -0.004 -0.045 
Fuel -0.05 0.108 1 0.151 0.193 0.178 0.113 
Subsistence -0. 15 0.372 0.151 1 0.213 -0.317 -0.09 
Ports & Misc 0.094  .0.204 0.193 0.213 1 0.173 0.639 
Equipage 0.388 -0.004 0.173 -0.317 0.173 1 0.249 
Voyage Repairs 0.396 -0.045 0.113 -0.09 0.639 0.129 1 


Figure 59 is the plot of salary versus voyage repairs. 
After the high leverage point is removed, the correlation falls to such a low 
level that it is considered insignificant. 

Figure 60 is the plot | of training versus subsistence cost 
after removal of the high leverage point. With the high leverage point 
removed, the correlation actually increases. In this case, a Spearman Rank 
Correlation of 0.894 will be used in the Cost Simulation Tool Model. During 
the cost analysis phase of the project, an explanation for this high correlation 
will be sought. 

Figure 61 is the plot of port and miscellaneous cost versus 
voyage repairs. With the high leverage point removed, the correlation 1s 


reduced to the point where it will be considered insignificant. 
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Figure 59 


T-ATF 166 Salary vs Voyage Repair Cost Scatter Plot 
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Figure 60 


T-ATF 166 Training vs Subsistence Cost Scatter Plot 
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Figure 61 
T-ATF 166 Port and Misc. vs Voyage Repair Cost Scatter Plot 
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V. THE MODEL 


The cost data was obtained from the Accounting Department at 
MSCPAC. This data was downloaded from the Financial Management 
Information System (FMIS). The format chosen for download was the 
Nucleus Report, a monthly report in which the costs incurred for each ship 
during the month are listed by cost category, with all sub-category line items 


listed. 
A. COST ANALYSIS MODEL 


The spreadsheet model is based on a normal costing system. The 
actual costs recorded monthly for the seven cost categories will be considered 
as direct costs since they are easily attributable to each ship using FMIS. 
Overhead costs are considered indirect costs, and budgeted costs will be used. 
Figure 62 1s an overview of the costing system used for the Cost Simulation 


Tool. 


Each run of the simulation will yield direct cost estimates for each of 
the seven cost categories used in the model. The simulation software 
generates a new value for each run as a result of a Monte Carlo simulation 
which will be discussed in detail in Chapter VI. The simulation generates a 
value based on the probability distribution assumed for each category. These 


assumptions are based on the data analysis. 
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Figure 62 
Cost Simulation Model Overview 
The seven direct cost categories are shown as triangles in Figure 62. 


These costs are represented as random variables for which a new value is 
produced with each replication of the simulation. 
summed by the spreadsheet. 


The values are then 


The two indirect cost categories are taken into account in the model as 


budgeted indirect overhead costs. The value is allocated as a pre-determined 


percentage of the sum of the direct costs for the replication. The indirect cost 
value 1s then summed by the spreadsheet. 
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B. MICROSOFT EXCEL SPREADSHEET MODEL 


The model was designed to capitalize on the Nucleus Report format. 
Figure 63 1s a printout of the T-AO 187 Class Cost Simulation Tool model. 
The model consists of a Microsoft Excel worksheet where the cost categories 
are listed and summed for a period of time. The user can choose the period of 
time desired and the number of ships to sum. The mean monthly amount for 
each cost category and a figure for the estimated amount for each cost 


category are listed for the assumed period of time and number of ships. 


The mean monthly amount for each category is based on the results of 
the data analysis covered in detail in Chapter IV. The user can change the 
mean amount using the edit assumptions button for the category. This 
allows the user to perform “what if’ analysis by increasing or decreasing the 
mean monthly amount used for the simulation run. When the Edit button is 
clicked, the Crystal Ball dialog box for the assumption appears. The user 
must be careful to enter his desired entry in the correct edit box. It is 
essential that the user review the Crystal Ball user's manual prior to using 
the spreadsheet model. The underlying distributions should not be changed 
by the user since these distributions were the result of extensive data 


analysis. 
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TAO Time Analysis Sheet 
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Figure 63 
TAO-187 Class Time Analysis Worksheet 
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$40, 340. 57 


|, Run | 


The user may also edit the assumed value for budgeted overhead costs. 
This value is not affected by Crystal Ball since it is deterministic. The user 
may change this value by increasing or decreasing the budgeted overhead by 


a percentage. 


The user may choose a time period to run the simulation for by 
clicking on the edit button for that assumption in the spreadsheet. The user 
must then enter a value for time m either days, weeks, or months. If the 


user does not select a time, the default value for time 1s one month. 


The user may choose the number of ships for the simulation run by 
clicking on the corresponding edit button. If the user does not provide an 


input, the default value for number of ships will be one. 


If it is desired to input a target operating income (profit), the user 
must click on the profit button at the bottom of the sheet. This will cause the 
total cost forecast to include the target income in its total. The target 
operating income can be in terms of a percentage of total cost or a specified 
amount depending on the desires of the user. The default setting for profit is 


zero since this 1s primarily a cost model. 


Once all of the options have been selected by the user, the simulation 
15 started by clicking the run button at the bottom of the sheet. This will 


cause the simulation to commence using whatever run preferences have been 
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input into Crystal Ball. If the user desires to change the run preferences, 
this is done using the pull down menu. The Crystal Ball user's guide 
contains detailed instructions for changing the dialog box that appears for 


run preferences. 


During the simulation run, a graphical representation of the empirical 
distribution created by the simulation values for total cost appears on the 
screen. Figure 64 1s a screen capture showing this graph. The graph allows 


the user to follow the progress of the simulation. 


At the completion of the simulation, similar graphs will be displayed 
for each of the cost categories. This allows the user to analyze each cost 
category separately. If the user desires a detailed report of the simulation 


results, he must click the “Report” button at the bottom of the spreadsheet. 
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VI. SIMULATION METHOD 


The purpose of this chapter is to present a brief description of the 
simulation method used by Crystal Ball, the simulation software used for 
this model. Crystal Ball uses a probabilistic Monte Carlo method to generate 


cost values which are summed on a Microsoft Excel spreadsheet. 


Problems handled by Monte Carlo methods are of two types, 
probabilistic or deterministic, depending on whether or not they are directly 
concerned with the behavior and outcome of random processes. In our case, 
the chief assumption has been that es direct costs behave according to 
probability distributions. Each cost category was fitted with a probability 


distribution in the data analysis performed in Chapter IV. 


In the probabilistic Monte Carlo case, the simplest approach is to 
observe random numbers, chosen in such a way that they directly simulate 
the random processes of the cost categories, and to infer the desired solution 
from the behavior of these random numbers. A probability density function 
is considered to have a total area, enclosed by the curve and the horizontal 
axis, of unity. A random number with a value between zero and one is 
generated by the computer. As seen in Figure 65, a point on the curve can 
be selected where the area beneath the curve accumulated to that point is 


equal to the value of the random number. 
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Атеа = х 
Random Number = x 


A = value generated by x 


Figure 65 
Random number relationship to distribution. 


There is a point, A on the x axis which is the projection of the point on 
the curve determined by the random number generation. The value of A 


corresponds to the value to be summed on the spreadsheet. 


This process is repeated for each of the cost categories using the 
probability distribution that was determined by the previous data analysis. 
After all values are obtained, the spreadsheet is summed to determine the 


total cost value for the period of time and number of ships concerned. 


The advantage of using a computer to perform these tasks is that the 
process can quickly be replicated to generate a new set of values. With each 
replication, another value is created for each cost category and the total cost. 
These values generated for each run are plotted so that a graph is generated 


which represents an empirical distribution of all the values generated for the 
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simulation. This empirical distribution is also analyzed to determine the 
mean and other salient characteristics of the values generated by the 
simulation. The graph can also be examined to find ranges of values which 


fit our desired certainty for the result. 
A. RANDOM NUMBER GENERATION 


An essential feature of Monte Carlo simulation is that at some point a 
substitution must be made for a random variable using a set of randomly 
generated values having the statistical characteristics of the random 
variable. The values that are ubi are called random numbers on the 
basis that they could well have been produced by chance by a suitable 
random process. As it turns out, the random numbers are not produced in 
this way, however, this should not affect our use of them. The question is not 
“Where did these numbers come from?” but rather “Are these numbers 
correctly distributed?” This question is answered by statistical tests on the 


random numbers themselves that are beyond the scope of this effort. 


When the term random number is used in Monte Carlo simulation, the 
standardized uniform distribution, U(0,1), is being referred to. In our case, 
the numbers are being generated by a pseudo-random number generation 
algorithm using a personal computer. The great advantage of this method is 
that the sequence of pseudorandom numbers can be exactly reproduced for 


purposes of computational checking. 
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Crystal Ball makes use of the Lehmer congruential method’. This 
method capitalizes on the architecture of the computer being used. In the 
case of the personal computers being used, all having 32 bit central 
processing units (CPU), the Lehmer method can generate a very long 
sequence of numbers without repeating a number. A 32 bit CPU permits 
4.25 bilhon numbers to be generated without a repetition. This application 


will never test these limits! 


B. TRANSFORMATION OF PSEUDO-RANDOM NUMBERS 


An important question 1s "How will the pseudo-random numbers from 
the uniform distribution be used to represent the various probability 
distributions in this model?” The method used to apply these pseudo-random 
variables to probability distributions is based on the inverse of the 
cumulative distribution function: 

Suppose that there are uniform U(0, 1) random variables 1), 15, ..., Un. 
If X, ~ F(x), we can generate X, by: 

X, 2 Fu) 


For F! to be unique, F must be strictly monotonic. 


The significance of this procedure is that random variables from 
known distributions can be represented by a sequence of uniform 
pseudo-random numbers if those numbers are "pushed" through the inverse 
‘Decisioneering Inc., Crystal Ball 3.0 Users Manual, p 38, Decisioneering 


Inc., Denver, CO. 
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cumulative distribution function of the known distribution being 


represented. 
C. CRYSTAL BALL 


Crystal Ball has sixteen distributions available for use. For each of 
these, Crystal Ball pushes a uniform pseudo-random number through the 
inverse transformation of the selected probability distribution. The resulting 
number is a random number from the selected distribution which is used for 


the present replication. 


Figure 66 is a simplified diagram that shows the sequence of events in 
the Crystal Ball simulation process. After the distributional assumptions 
are made and the values and parameters are entered for the spreadsheet 
cells in which the simulation values are generated (referred to by Crystal 
Ball as the assumption cell), the simulation is started. Next, the random 
numbers are generated for each assumption cell and the transformations are 


made. 


Once the simulated values for each assumption cell are obtained, the 
spreadsheet is calculated. The result of this calculation yields the forecast 
values being sought. Once the forecast values are obtained, the process is 
repeated until the number of runs entered by the user is reached. On each 
run, the values for each category of cost are averaged with all of the previous 


values obtained in the previous runs. This average, or grand mean, is the 
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Simulation Cycle 


Generate 
Assumption Values 


Recalculate 
Spreadsheet 


Extract 
Forecast Results 


9 


Figure 66 


Crystal Ball simulation overview. 


mean value achieved so far for the simulation. At the end of the simulation, 
the mean value represents the final grand mean value for the entire 


simulation. 
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VII. VALIDATION OF DATA ANALYSIS 

The distributional assumptions made in Chapter IV were based on 
using a partition of the entire data set consisting of two chronological halves. 
The earlier half of the entire monthly cost data set was used during the data 
analysis. The later half of the entire data set was next used to validate the 
assumptions made in Chapter IV. If the distributional assumptions made 
earlier remain valid, the validation data set should be well fit by the same 
distributions as those chosen earlier. 

An additional check was made by running the Cost Simulation Tool 
model against several cases of the actual data for randomly selected months. 
If the data analysis and model are valid, the simulation run values should 
reasonably represent the actual monthly values. 

A. SUMMARY OF DATA ANALYSIS VALIDATION RESULTS 

The validation data sets were analyzed in two ways. First, probability 
distribution fit plots were made for each cost category to assess the fit of the 
distribution family found for the corresponding cost category in Chapter IV to 
the validation data. For these fits, the parameters were estimated using the 
validation data. A second probability distribution fit was also attempted 
using the same distribution family and the parameters estimated in Chapter 


IV. 
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For every cost category, a good fit to the validation data was obtained 
with the distribution family used for the corresponding cost category in 
Chapter IV, but with new parameter estimates. When attempting to fit the 
validation data with the parameters (using same distributions) found in 
Chapter IV, the fits obtained were almost always poor. This result caused 
another look to be taken at the data. Since the sample size was relatively 
small with only half of the data set, both halves were used to get more precise 
estimates of the parameters. 

This time, the entire data set (filtered as in Chapter IV) for each 
category was used. Again, two fits were made for each cost category. The 
first set of plots was made by attempting to fit the data set with both the 
distribution and parameters from Chapter IV. The second set of plots was 
made by attempting to fit the data set with only the distribution from 
Chapter IV, allowing AGSS to automatically use parameters that provide the 
best fit . In this round of plotting, the fits to the full data set were better 
than those found with the validation data set when using both the 
distribution and parameters from Chapter IV as would be expected. 
Excellent results were achieved when fitting the distribution only to the full 
data sets. The model was updated to reflect the new results obtained from 


allowing AGSS to automatically select parameters that best fit the entire 
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data set. Validation test runs of the model were again made after the 
changes were made to the model. 

Figure 67 shows the original probability distribution plot for T-ATF 
166 Class Voyage Repairs which was based on the first half of the original 
partitioned data set. Figure 68 is the new probability distribution plot for 
the same cost category using the full data set. It is evident from comparing 
the two views that the data is equally well fit by the same distribution 


though the parameters have changed slightly. 


T-ATF 166 Class Voyage Repair Cost 


LOGNORMAL DENSITY FUNCTION, N=16 LOGNORMAL CUMULATIVE DISTRIBUTION FUNCTION, N=16 
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0 210* 410* 610% о 210* 410* 610* 
Dollars Dollars 


LOGNORMAL PROBABILITY PLOT, N=16 





Figure 67 
T-ATF 166 Class Voyage Repair Probability Distribution Fit Plot. 
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T—-ATF 166 Voyage Repair Cost Validation Data 
LOGNORMAL DENSITY FUNCTION, N-29 LOGNORMAL CUMULATIVE DISTRIBUTION FUNCTION, N=29 
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Figure 68 
T-ATF 166 Revised Voyage Repair Probability Distribution Fit 


All cost categories were treated in the same way and the results were 
compared to insure that the two sets of parameters found (Chapter IV and 
Chapter VII) for each cost category were within a 95 percent confidence 
interval of each other (based on fitting 95 percent confidence interval to 
Chapter IV estimates and testing to ensure that Chapter VII estimates fell 
inside.) All data sets passed this test. A brief summary of the results of the 
improved distribution fits obtained using the entire data sets is given in 


Tables 5 and 6. 


Table 5 


T-AO 187 CLASS DATA ANALYSIS VALIDATION RESULTS 


Cost Category Distribution 


Salary 
Training 

Fuel 
Subsistence 
Port and Misc. 
Ship's Equipage 
Voyage Kepairs 


Logistic 
Weibull 
Logistic 
Logistic 
Weibull 
Weibull 
Gamma 


Parameters 
o=732,000, B293,164 
C=1.5962, a=19,615 
о=20$ А64, B=48,002 
01927307 5—31570 
C- 1.579, a-220,080 
C-0.96836, a- 19,117 
0-097569, B2119,370 


Table 6 


Mean StdDev Fit Quality 


125,644 
17,689 
204,943 
19,342 
196,010 
19,492 
118,400 


170,648 
12,008 
105,861 
1,158 
141,300 
21,107 
971767 


Good 
Good 
Good 
Fair 
Fair 
Good 
Good 


T-ATF 166 CLASS DATA ANALYSIS VALIDATION RESULTS 


Cost Category Distribution 
Normal 
Gamma 
Weibull 
Logistic 

Lognormal 

Lognormal 

Lognormal 


Salary 

Training 

Fuel 
Subsistence 
Port and Misc. 
Ship's Equipage 
Voyage Repairs 


Parameters 


и=143,830,о=74.280 
a-1.791, 821,44 
C=1.2947, a=60,386 56,961 
a=3,424, B=9162 
u=10371,0=1.0614 
u=75506, o=0.4695 
и=11236,о=0.95045 


2,588 


3,642 
59,353 
4,111 


115,380 


14.280 
2,356 
49,056 
1,902 
83,532 
5,274 
108,540 


Mean Std Dev Fit Quality 
143,830 


Fair 
Good 
Fair 
Marginal 
Good 
Good 
Good 


B. VALIDATION RESULTS USING COST SIMULATION MODEL 


Once the model parameters were updated as a result of the previous 
graphical data analysis, the Cost Simulation Tool was run for several cases 
to compare the results of the simulation to the actual values documented 


with historical cost data. These validation test runs will be discussed on an 


individual case basis. 
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1. Case 1: T-AO 187 Class Per Diem Rate Comparison 

MSCPAC charges their customers for services provided by their 
ships on a per diem basis. Each ship class has a per diem rate that has been 
computed by MSC headquarters and is charged to the sponsor for each ship 
day of use (this per diem rate includes overhead costs and the corresponding 
model used for the validation runs also included overhead costs). In this 
case, the T-AO 187 Class model was run for a time period of one ship day. 
The results were compared against the per diem rate of $60,765 per ship day 
(FY 94 rate). If the model and data analysis accurately represent the actual 
costs incurred, the result obtained by running the model for one ship day 
should be close to the per-diem rate presently being charged. 

Figure 69 15 а combined empirical — probability 
distribution/frequency plot generated by the simulation. The actual rate of 
$60,765 falls at the 98th percentile of the empirical distribution. This 
implies that approximately 98 percent of the time, the actual per diem rate 
charged will be higher than the result of a simulation run. This result 
indicates that either the model is predicting lower values than actually 
necessary to cover operating costs, or the present per diem rate is set too 
high. Since the per diem rate falls within the values generated by the 


simulation, the simulation method does seem to generate values similar in 
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Forecast: Total Cost Forecast 









Frequency Chart 1,989 Trials Shown 
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Certainty Range is from -Infinity to $60,765.00 Dollars 





Figure 69 
T-AO 187 Class Per Diem Comparison Result. 


value to the per diem rate. The next case, which simulates T-AO 187 direct 
costs only, provides some insight. 
2. Case 2: T-AO 187 Class Historical Data Comparison Method 

In this case, the Cost Simulation Tool was run for a period of 
one month for one ship. The results of this simulation run were then 
compared to the actual monthly costs reported for all ships (data set was 
filtered according to the model assumptions). Table 7 shows a randomly 
selected actual monthly data set for five ships and indicates how many 
standard deviations away from the simulation mean for the cost category 
(from Table 5) the selected data was. A better view of the model accuracy is 
seen in Figure 70 that compares the actual total costs with the simulation 


mean total cost. An interval of one standard deviation of the total cost 
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Table 7 
T-AO 187 CLASS COMPARISON (ACTUAL VS SIMULATION) 


Ship 1 #50 Ship 2 #50 Ship 3 #50 Ship 4 #50 Ship 5 #50 
Cost Category 


Salary 319,177 -2.382 510,941 -1.299 793,971 0.383 961,026 1.370 964,076 1.388 
Training 11,081 -0.550 17,937 0.056 3.251 -1.230 51.928 2976 22,291 0.437 
Fuel 243,552 0.365 415,734 2.610 169,215 -0.564 126,511  -1.114 140,132  -0.938 
Subsistence 29,651 1.647 18,699 -0.110 29,651 1.647 20,461 0.173 11,948 -1.193 
Port and Misc. 586,889 2.766 115,181 -1.016 196,304 0.046 280,271 1.145 186,436 -0.084 
Ship's Equipage 43,601 1.082 154 -0.896 16,213 -0.183 5,559 -0.675 13,086 -0.327 
Voyage Repairs 264,090 1.490 312,193 1.565 202048 0667 425,969 2.492 113,970 -0.051 
Total Cost 1,498,041 0.813 1,390,739 0.349  1.410,653 0.435 1,871,725 2.429 1,451,939 0.614 

2500 
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Data Points 
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Figure 70 
T-AO 187 Class Monthly Actual Data Comparison Result. 


obtained by the model run above and below the simulation mean is also 
shown 1n the plot. 

For the monthly data comparison run, the budgeted overhead 
feature of the model was turned off. This was done because there is no actual 
overhead data included with the actual cost data as reported by FMIS. The 
overhead feature is used in the model for estimating the per-diem rate only; 


when comparing to actual cost categories listed in FMIS, there is no category 


109 


for budgeted overhead. The actual costs are mostly (69.496) within a range of 
one standard deviation of the simulation mean value. 

The observation that the model seemed to predict values lower 
than the actual per-diem rate (as seen in case 1) would imply that the model 
would predict lower monthly totals than the actual monthly data. This is 
not seen in case 2. The model values seem to be about the same as the actual 
values. The actual direct costs are all close to the simulation mean value 
(with the budgeted overhead feature turned off). This suggests that either 
the budgeted overhead rate used in the model is lower than what would be 
needed to make the actual per-diem rate accurate or that the per-diem rate is 


too high . It appears that the model is predicting direct costs accurately. 


3. Case 3: T-ATF 166 Class Per Diem Rate Comparison 
In this case, the per diem rate being charged for FY 94 is 
$15,955. Figure 71 is an empirical probability distribution/frequency plot 
generated by the simulation. This plot indicates that the actual per diem 
rate falls at the 77.6 percentile. This result indicates that the model 
predictions are more accurate than those of case one and that the actual per 
diem rate being charged for this class vessel is closer to the actual costs of 


operation than was found in case one. 
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Figure 71 
T-ATF 166 Class Per Diem Comparison Result. 
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4. Case 4: T-ATF 166 Class Historical Data Comparison Method 


This is similar to case two discussed in section B above. Table 8 


is the actual monthly cost comparison data by cost category for five randomly 


selected T-ATF 166 class ship-months. 


Other than a few isolated instances, the actual data was close to 


the simulation mean values in every cost category (from Table 6). In most 


Cost Category 
Salary 

Training 

Fuel 
Subsistence 
Port and Misc. 
Ship's Equipage 
Voyage Repairs 
Total Cost 


Ship 1 


109,516 
12,178 
43,866 

3, 459 
9,289 
4, 159 
132,368 
314,835 


#50 


-0.535 
9.157 
-0.283 
-0.078 
-0.579 
0.030 
0.156 
-0.464 


Ship 2 


130,546 
1,783 
44,584 
5,771 
7,674 
25,619 
122,597 
338,574 


Table 8 
T-ATF 166 CLASS COMPARISON (ACTUAL VS SIMULATION) 


#SD 


-0.024 
-0.407 
-0.267 
1.406 
-0.598 
4.527 
0.067 
-0.304 


Ship 3 


131,804 
481 
11,671 
219 
94,103 
10,486 
154,183 
405,647 
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#50 


-0.224 
-1.104 
-1.004 
-0.424 
0.409 
1.354 
0.355 
0.149 


Ship 4 


141,973 
1,756 
10,042 
3,441 
14,431 
860 
23,782 
256,285 


"USD 


-0.081 
-0.422 

0.303 
-0.089 
-0.519 


-0.662 
-0.835 
-0.860 


Ship 5 


259,192 
3,519 
27,509 
1,361 
32,801 
9,191 
285,624 
619,197 


#50 


1.559 
0.522 
-0.650 
-1.424 
-0.305 
1.084 
1.554 
1.590 


cases, the values were within one standard deviation of the simulation mean 
value. Figure 72 is a plot that compares the actual total costs with the 
simulation mean total cost. The values are mostly (72.4%) within one 
standard deviation of the simulation mean. The simulation produces a 
reasonable representation of the data from this ship class . 
5. Summary of Model Comparison Results 

The results seen in cases three and four suggest that the model 
concept is correct but that the overhead assumptions in the T-AO 187 class 
model need to be examined for accuracy. In both cases, the majority of the 
actual results were within one standard deviation of the simulation mean. 
As the primary focus has been the simulation of actual values, the model is 


viable in that respect. 
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Figure 72 
T-ATF 166 Class Monthly Actual Data Comparison Result. 


112 


VIII. CONCLUSIONS 

The original intent of this project was to develop a tool that would 
enable MSCPAC to perform reasonably accurate "what if" cost analyses for 
the ships they own and operate. The Cost Simulation Tool was developed to 
fill that requirement. As in any computer application, the expected results 
from the product are only as good as the inputs to the program. This is the 
famous "garbage in = garbage out" principle of computing. 

A qualified success has been achieved with the Cost Simulation Tool. 
The success is in the ability of the simulation to accurately forecast direct 
ship operating costs. Unfortunately, in the area of indirect overhead costs, 
more work is needed, which is beyond the scope of this thesis. The overhead 
estimate used in the Cost Simulation Tool model was a best guess figure 
provided by MSCPAC operations department and is not documented 
anywhere. To correct this problem would require a cost analysis to 
determine the extent of the overhead costs for MSC. The total overhead costs 
might then be allocated as a percentage of total direct costs. 

The original problem statement expressed by the MSCPAC comptroller 
included some degree of certainty about the costs in the areas of 
infrastructure, overhead, and chartering services. As this project has 
unfolded, the concerns expressed by the comptroller in the area of direct costs 


for the ships owned and operated by MSC were real, but not as serious as 
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perceived. This analysis has shown that the cost information available from 
the FMIS was accurate but cumbersome to access and use. The Cost 
Simulation Tool has been demonstrated to quickly and reliably produce 
operating cost estimates for the direct costs. The overhead costs; however, 
are not precisely known, and when the rough estimate is included in the cost 
simulation model, the accuracy of the model suffers. 

The problem of imprecise overhead estimates is not unique to MSC. 
Since the advent of DBOF, most Department of Defense activities subjected 
to this mode of funding have suffered from an inability to accurately forecast 
overhead amounts to be included in billing. The subject of overhead cost 
analysis is wide open and needs further examination in the case of MSC. 

The direct costs of operating the ships are reported in great detail and 
can be directly traced to their source in an economically feasible manner. 
The Financial Management Information System (FMIS) maintains records of 
every voucher written by the ships. Clearly, the record keeping is complete 
in this area. The Cost Simulation Tool exploits this historical database to 
produce accurate direct cost results by simulation. 

A project of enormous use to MSC management would be to enhance 
FMIS to include a query system that could allow the managers to access data 
and manipulate it from their Microsoft Excel spreadsheets. This would 


permit the data of interest to be displayed in tables and graphs without the 
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necessity of mastering the complexities of the FMIS system itself. The 
combination of the Cost Simulation Tool and the query system would give the 
managers great flexibility in manipulating the data and analyzing the 
outcome of various scenarios. The two systems could work together to give 
management a better understanding of the operating costs of their business. 

Another analysis could be conducted to determine the overhead costs 
and their allocation at MSC. The result of this future analysis could be 
incorporated into the Cost Simulation Tool to improve its accuracy. The Cost 
Simulation Tool could then be used to calculate per diem rates that would 
reflect actual MSC costs with much greater accuracy. 

Without following up to determine the overhead costs and their 
drivers, use of the Cost Simulation Tool is limited to estimating direct costs 
only. The Cost Simulation Tool was shown to produce accurate results in the 
area of direct costs, but the total cost accuracy is subject to the rough 
estimate of budgeted overhead. This is because the overhead cost calculation 
is, at present, based on a "best guess" approximation of overhead costs. 

The question is, "What has been accomplished by the Cost Simulation 
Tool?" First, the data analysis has shown that the operating cost data is 
indeed random in nature and can be represented by probability distributions. 
Second, once the probability distributions have been determined, the costs 


can be simulated using the Monte Carlo simulation method. 
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These results, when combined, give the cost accountant a new tool to 
use in forecasting the costs of doing business. No longer tied to regression 
alone; through simulation, the cost accountant can see the entire range of 
behavior for the costs of concern. This ability can give forecasts much greater 


accuracy than in the past. 
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APPENDIX A 


Figure 2 Data Table 


ANALYSIS OF LOGISTIC OISTRIBUTION FIT 


OATA OSALRY14 
SELECT ION ALL 
X AXIS LABEL: Оо1ісгз х 964,075 
SAMPLE SIZE : 25 
CENSOR ING : NONE 
FREQUENCIES : 1 
EST. METHOD - MAXIMUM L KELIHOOO 
CONF METHOD : ASYMPTOTIC NORMAL APPROXIMATION 
CONF. INTERVALS COVARIANCE MATRIX OF 
(95 PERCENT) PARAMETER ESTIMATES 
PARAMETER ESTIMATE LOWER UPPER ALPHA BETA 
ALPHA 0.71841 0.65081 0.78602 0.0011892 0 
BETA 0.099551 0.06691 0.13219 0 0.00027722 
LOG LIKELIHOOD FUNCTION AT MLE = 8.0678 
SAMPLE FITTED GOODNESS OF FIT TESTS 
MEAN 0.70757 0.71841 CHI-SQUARE : 2.2654 
STO OEV : 0.17471 0.180$7 OEG FREED: 1 
SKEWNESS:  —0.34818 0 SIGNIF 2 021222 
KURTOSIS: 2.2898 4.2 KOLM-SMIRN : 0.11668 
« BASED ON MIOPOINTS OF FINITE INTERVALS SIGNIF 0.88546 
CRAMER-V M 0.0632SS 

PERCENTILES SAMPLE FITTED SIGNIF > .15 

5: 0.40366 0.42529 ANDER-DARL : 0.4681 

10: 0.4413 0.49968 SIGNIF > 215 

25: 0.61702 0.60904 

50: О. 75042 0.71841 KS, AD, AND CV SIGNIF. LEVELS NOT 

7S: 0.823S6 0.82778 EXACT WITH ESTIMATED PARAMETERS. 

90: 0.89774 0.93715 

95: 0.99684 1.0115 NOTE: A SMALL SIGNIFICANCE LEVEL 


CHI-SQUARE GOODNESS OF FIT TABLE 


(EG. Ps.01) INDICATES LACK OF FIT 


LOWER UPPER OBS EXP 0-Е ((0-Е).2)=Е 
-INF. 0.61714 7 6.6388 0.36123 0.019655 
0.61714 0.71999 3  5.960S -2.960S 1.4705 
0.71999 0.82285 8 5.9157 2.0843 0. 73436 
0.82285 +INF. 7 6.485 0.515 0.040899 
TOTAL 25 25 2.2654 
Figure 7 Data Table 
ANALYSIS OF WEIBULL OISTRIBUTION FIT 
OATA OTRA I NO7 
SELECTION : ALL 
X AXIS LABEL: Dollars 
SAMPLE SIZE : 2S 
CENSORING |: NONE 
FREQUENCIES : 1 
EST. METHOD : MAXIMUM L!KELIHOOO 
CONF METHOO : ASYMPTOTIC NORMAL APPROXIMATION 
CONF. INTERVALS COVARIANCE MATRIX OF 
(9S PERCENT) PARAMETER ESTIMATES 
PARAMETER ESTIMATE LOWER UPPER C a 
C (SHAPE) 1.5615 1.0715 2.0515 0.062478 2.5573Е2 
a (SCALE) 24854 18307 31422 258770 1.1189Е7 


LOC LIKELIHOOD FUNCTION AT MLE = –272.27 


SAMPLE FITTED GOODNESS OF FIT TESTS 
MEAN 22411 22346 CHI-SQUARE : 1.2157 
STD OEV : 14578 14622 OEG FREED: 2 
SKEWESS : 0.61912 1.0009 SIGNIF 0.54454 
KURTOSIS: 2.364 1.4288 KOLM-SMIRN : 0.09619 
* BASED ON MIOPOINTS OF FINITE INTERVALS SIGNIF : 0.97484 
CRAMER-V M : 0.031256 
PERCENTILES SAMPLE FITTEO SIGNIF 5 > -15 
S: 3281.4 3710.8 ANDER-DARL : 0.24434 
10: $637 5884.1 SIGNIF > 15 
25: 11395 11196 
50: 17937 19663 KS, AD, AND CV SIGNIF. LEVELS NOT 
75: 27069 30650 EXACT WITH ESTIMATED PARAMETERS. 
90: 45251 42418 
9S: 49853 50204 NOTE: A SMALL SIGNIFICANCE LEVEL 


(ЕС. Р<.01) INDICATES LACK OF FIT 


CHI-SQUARE GOOONESS OF FIT TABLE 


LOWER UPPER OBS EXP 0-Е ((0—).2)-Е 

-INF. 8473.4 3 4.2473 -1.2473 0.36631 

8473.4 16947 6 6.3229 —0.32287 0.016487 
16947 25420 7 5.5502 1.4498 О. 37869 
25420 33894 3 3.9125 -0.94252 0.22533 
33894 *INF. 6 4.937 1.063 0.22886 
TOTAL 25: 25 1.2157 
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Figure 11 Data Table 


ANALYSIS OF LOGISTIC DISTRIBUTION FIT 


OATA OFUEL 14 

SELECT I ON ЕГА ЕЕ 

X AXIS LABEL: Oollars X 534,166 
SAMPLE SIZE : 25 

CENSOR I NG : NONE 


FREQUENCIES 251 
EST. METHOO : MAXIMUM LIKEL!400D 
CONF METHOO : ASYMPTOTIC NORMAL APPROX IMATION 


CONF. INTERVALS COVARIANCE MATRIX OF 


(95 PERCENT) PARAMETER ESTIMATES 
PARAMETER ESTIMATE LOWER UPPER ALPHA BETA 
ALPHA 0.37821 0.30979 0.44664 0.0012183 0 
BETA 0.10076 0.067723 0.1338 0 0.000284 
LOG LIKELIHOOD FUNCTION AT MLE - 6.0708 
SAMPLE FITTED GOOONESS OF FIT TESTS 
MEAN : 0.38495 0.37821 KOLM-SMIRN : 0.15724 
STO OEV : 0.20786 0.18276 SIGNIF : 0.56679 . 
SKEWNESS: 0.72848 0 СКАМЕК-У М : 0.10111 
KURTOSIS: 5.1271 4.2 SIGNIF : >.15 
* BASEO ON MIDPOINTS OF FINITE INTERVALS ANDER-DARL : 0.71454 
SIGNIF : а 
PERCENTILES SAMPLE FITTEO 
5: 0.0063472 0.081531 KS, AO, ANO CV SIGNIF. LEVELS NOT 
10: 0.072774 0.15682 EXACT WITH ESTIMATED PARAMETERS. 
25 0.31678 0.26752 : 
50: 0.38108 0.37821 NOTE: A SMALL SIGNIFICANCE LEVEL 
753 0.46914 0.48891 (EG. Ps.01) INDICATES LACK OF FIT 
90: 0.52469 0.59961 
95: 0.77829 0.6749 
Figure 15 Data Table 
ANALYSIS OF LOGISTIC DISTRIBUTION FIT 
OATA : OSUBSO9 
SELECTION : ALL 
X AXIS LABEL: Dollars X 40177 
SAMPLE SIZE : 25 
CENSORING : NONE 


FREQUENCIES : 1 
EST. METHOO : MAXIMUM LIKELIHOOD 
CONF METHOD : ASYMPTOTIC NORMAL APPROXIMATION 


CONF. INTERVALS COVARIANCE MATRIX OF 


(95 PERCENT) PARAMETER ESTIMATES 
PARAMETER ESTIMATE LOWER UPPER ALPHA BETA 
ALPHA 19486 16957 22016 1.6651Е6 0.0000Е0 
ВЕТА 372591 2503.7 4946.4 0.0000Е0 5.8815Е5 
LOG LIKELIHOOO FUNCTION AT MLE - -257.48 
SAMPLE FITTED GOOONESS OF FIT TESTS 
MEAN : 19824 19486 KOLM-SMIRN : 0.21565 
STO OEV : 7997.5 6756.5 SIGNIF : 0.19532 
SKEWNESS: 0.44477 О CRAMER-V M : 0.27207 
VURTOSIS: 5.0097 4.2 SIGNIFSEES — 15 
* BASED ON MIDPOINTS OF FINITE INTERVALS ANDER-DARL : 1.5124 
SIGNIF : > .15 
PERCENTILES SAMPLE FITTEO 
S: 5870.5 8518 KS, AO, ANO CV SIGNIF. LEVELS NOT 
10: 13177 11301 EXACT WITH ESTIMATED PARAMETERS. 
25: 17996 15394 
БО: 19749 19486 NOTE: A SMALL SIGNIFICANCE LEVEL 
75: 20461 23579 (EG. Р<.01) INDICATES LACK OF FIT 
90: 29651 27671 
95: 38236 30454 
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Figure 19 Data Table 


ANALYSIS OF WEIBULL OISTRIBUTION FIT 


OATA OPORTOS 

SELECT I ON ALL 

X AXIS LABEL: Dollors 

SAMPLE SIZE : 25 

CENSOR I NG : NONE 

FREQUENCIES : 1 

EST. METHODO : MAXIMUM LIKELIHOOD 

CONF METHOO : ASYMPTOTIC NORMAL APPROXIMAT ION 


PARAMETER ESTIMATE LOWER 


CONF. INTERVALS COVARIANCE MATRIX OF 
(95 PERCENT) PARAMETER ESTIMATES 
UPPER C 


a 
0.18694 1.7119ЕЗ 


C (SHAPE) 2.7352ЕО 1.8876Е0 3.5828ЕО 
a (SCALE) 1.7101ES 1.4530ES 1.9673E5 1711.9 1.7205E8 
LOG LIKELIHOOD FUNCTION AT MLE = —310.59 
SAMPLE FITTED GOOONESS OF FIT TESTS 
MEAN 1.5257E5 1.5215Е5 CHI-SQUARE : 0.4684 
STO OEV : 6.1177E4 6.0071Е4 OEG FREEO: 2 
SKEWNESS: 8.4080Е-2 6.3949Е-1 SIGNIF 0.7912 
KURTOSIS: 2.9129Е0 1.2909Е0 KOLM-SMIRN : 0.089326 
e BASED ON MIDPOINTS OF FINITE INTERVALS SIGNIF : 0.98844 
CRAMER-V M : 0.026164 
PERCENTILES SAMPLE FITTEO SIGNIF END SS 
Sc 6.2616Е4  5.7733E4 ANDER-DARL : 0.1888 
10: 7.4355E4  7.5114E4 SIGNIF > .15 
25: 1.1518Е5 1.0845Е5 
50: 1.5866E5 1. 4957Е5 KS, AD, AND CV SIGNIF. LEVELS NOT 
292 129297655  12927ЛЕӘ EXACT WITH ESTIMATEO PARAMETERS. 
90: 2.2006ES . 2.3199E5 
95: 2.6501Е5 2.5541Е5 NOTE: A SMALL SIGNIFICANCE LEVEL 


(EG. P<.01) INDICATES LACK OF FIT 


CHI-SQUARE GOOONESS OF FIT TABLE 


LOWER — UPPER OBS EXP O-E . ((0-£)«2)«€ 
—INF. 8.7443E4 4 3.6893 0.31067 0.026161 
8.744564 1.3116Е5 5 5.9051 -0.9051 0.15816 
1.3116Е5 1.74895 8 6.7734 1.2267 0.22214 
1.7489Е5 2.1861Е5 5 5.1033 -0.10326 0.0020894 
2.1861E5 *INF. > 3.531 -0.53096 0.079841 
TOTAL 25 25 0.4684 
Figure 23 Data Table 
ANALYSIS OF WEIBULL DISTRIBUTION FIT 
OATA OE£QUIPOS 
SELECT I ON ALL 
X AXIS LABEL: Dollars 
SAMPLE SIZE : 25 
CENSOR I NG : МОМЕ 
FREQUENCIES : 1 
EST. METHOD : | MAXIMUM LIKELIHOOD 
CONF METHOD : ASYMPTOTIC NORMAL APPROXIMATION 
CONF. INTERVALS COVARIANCE MATRIX OF 
(95 PERCENT) PARAMETER ESTIMATES 
PARAMETER ESTIMATE LOWER UPPER C a 
C (SHAPE) 1.0816 О. 73862 1.4247 0.030616 2.1536E2 
a (SCALE) 21037 13039 29035 215.36 1.6646Е7 
LOG LIKELIHOOD FUNCTION AT MLE = -273.04 
SAMPLE FITTED GOOONESS OF FIT TESTS 
MEAN 20456 20414 CHI-SQUARE : 1.4944 
STD DEV : 17708 18890 DEG FREED: 1 
SKEWNESS: 1:0125 152119 SIGNIF : 0.22154 
KURTOSIS: 3.1397 1.6612 KOLM-SMIRN : 0.10514 
* BASED ON MIDPOINTS OF FINITE INTERVALS SIGNIF : 0.94511 
СКАМЕН-У М : 0.044193 
PERCENTILES SAMPLE FITTED SIGNIF с 2:015 
Si 755-5 1350.2 ANDER-DARL : 0.29121 
10: 1958 2626.8 SIGNIF >: 15 
25: 9331.6 6648.7 
50: 14067 14991 KS, AD. AND CV SIGNIF. LEVELS NOT 
75: 28376 28453 EXACT WITH ESTIMATEO PARAMETERS. 
90: 46797 45484 
93: 50740 58012 NOTE: A SMALL SIGNIFICANCE LEVEL 


LOWER 
-INF. 
11001 
22001 
33002 
TOTAL 


(EG. P<.01) INDICATES LACK OF FIT 


CHI—-SQUARE GOODNESS OF FIT TABLE 


UPPER 085 ЕХР Ge 

11001 B 9.7754 17754 

22001 3 6.4733 2.5267 

33002 3 5.8408 -0.84081 

*INF. 5 4.9105 0.089504 
25 25 


((0-Е)»2)-Е 
0.52244 
0.98622 
0.18407 
0.0016314 
1.4944 


Tug 
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Figure 27 Data Table 


ANALYSIS OF GAMMA OISTRIBUTION FIT 


DATA ОУКЕРО6 
SELECT I ON ALL 
X AXIS LABEL: Ocllors 
SAMPLE SIZE : 25 
CENSOR | NG : NONE 
FREQUENCIES : 1 
EST. METHOD : MAXIMUM LIKELIHOOD 
CONF METHOD : ASYMPTOTIC NORMAL APPROXIMATION 
CONF. INTERVALS COVARIANCE MATRIX OF 
(95 PERCENT) PARAMETER ESTIMATES 
PARAMETER ESTIMATE LOWER UPPER ALPHA BETA 
ALPHA 8.3755E-1 0.43494 1.2402E0 0.042177 —6.4306E3 
BETA 1.2770Е5 45473 2.0993Е5 -6430.6 1. 7593Е9 


LOG LIKELIHOOD FUNCTION AT MLE = -339.34 


SAMPLE FITTED GOODNESS OF FIT TESTS 
MEAN 1.0695E5 1.0695Е5 CHI-SQUARE : 0.76465 
STO OEV : 9.7542Е4 1.1687Е5 OEG FREED: 1 
SKEWNESS: 9.3881Е-1 2.1854E0 SIGNIF 0.38188 
KURTOSIS:  3.4970E0 1.0164Е1 KOLM-SMIRN : 0. 11317 
* BASED ON MIDPOINTS OF FINITE INTERVALS SIGNIF 0.90603 
CRAMER-V M : 0.072187 
PERCENTILES SAMPLE FITTEO SIGNIF > .15 
5: 3.3061E3 3.3838ЕЗ ANDER-DARL : 0.47451 
10: 6.1625E3 7.8749Е3 SIGNIF 2.15 
25: 2.1116Е4  2.5245E4 
50: 8.3043Е4 6.8580Е4 KS, AD, AND CV SIGNIF. LEVELS NOT 
75: 1.8079Е5 1.4805Е5 EXACT WITH ESTIMATED PARAMETERS. 
90: 2.2253E5 2.5720Е5 
95: 2.2788Е5 3.4129Е5 NOTE: A SMALL SIGNIFICANCE LEVEL 


(EG. P<.0°) INOICATES LACK OF FIT 
CHI-SQUARE GOOONESS OF FIT TABLE 


LOWER 
UNES 


6.3591E4 


1.2718 
1.9077 
TOTAL 


DATA 

SELECT |I ON : 
X AXIS LABEL: 
SAMPLE SIZE : 
CENSOR I NG : 
FREQUENCIES : 
ЕЗІ МЕРНО : 
CONF METHOD : 


UPPER 
6.3591Е4 
1.2718Е5 
1.9077Е5 
TINF. 


ES 
ES 


OBS 
11 


EXP 


0-Е 


11.924 -0.9241 


5 5.6027 -0.6027 
2 11205 -0212053 


6 4.3526 


25 25 


1.6474 


((О-Е)»2)+Е 


0.071616 
0.064833 
0.0046631 
0.62354 
0.76465 


Figure 31 Data Table 


ANALYSIS OF NORMAL DISTRIBUTION FIT 


TSALRYO2 
ALL 
Dollars 

16 

NONE 

1 

MAXIMUM LI 
EXACT 


CONF. 
(95 


KEL I HOOO 


INTERVALS 


PERCENT) 


PARAMETER ESTIMATE LOWER UPPER 


MU 
SIGMA 


1.2388Е5 93687 
5.4864E4 41857 8.7699Е4 


1.5408Е5 


COVARIANCE MATRIX OF 


PARAMETER ESTIMATES 


MU 


SIGMA 
1.8813E8 0.0000ЕО 
0.0000Е0 9.4066Е7 


LOG LIKELIHOOD FUNCTION AT MLE = -197.31 


SAMPLE FITTEO GOOONESS OF FIT TESTS 
MEAN 1.2388E5 1.2388E5 KOLM-SMIRN : 0.15001 
STO OEV : 5.6664Е4 5.4864E4 SIGNIF : 0.86422 
SKEWNESS:  1.1985E-1 0.0000Е0 CRAMER-V M : 0.075901 
KURTOSIS: 4.1197ЕО 3.0000E0 SIGNIF : > 18 : 
» BASED ON MIDPOINTS OF FINITE INTERVALS ANDER-DARL : 0.47923 
SIGNIF > .15 
PERCENTILES SAMPLE FITTEO 
5: 1.3497E3 3.362054 KS, AO, ANO CV SIGNIF. LEVELS NOT 
10: 6.2848E4 5.356364 EXACT WITH ESTIMATED PARAMETERS. 
25: 9.396984 8.6895Е4 
50: 1.339365 1.2388E5 NOTE: A SMALL SIGNIFICANCE LEVEL 
75: 1.4800E5 1.6087Е5 (EG. Р<.01) ІМОІСАТЕ5 LACK OF FIT 
90: 1.6765Е5 1.9421Е5 
95: 2.5919Е5 2.1415Е5 
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Figure 35 Data Table 


ANALYSIS OF GAMMA OISTRIBUTION FIT 


OATA TTRAINO7 

SELECT I ON З АЦЕ 

X AXIS LABEL: Oollars 

SAMPLE SIZE : 16 

CENSOR ING : NONE 

FREQUENCIES : 7 

EST. METHOD : MAXIMUM LIKELIHOOO 

CONF METHOO : ASYMPTOTIC NORMAL APPROX! MATION 


CONF. INTERVALS COVARIANCE MATRIX OF 
(95 PERCENT) PARAMETER ESTIMATES 
PARAMETER ESTIMATE LOWER UPPER ALPHA BETA 
ALPHA 1.5843: 081501 222537 0.19665 —35.1855E2 
BETA 2242.5 552.64 3932.4 =318.55 7.4307E5 
LOG LIKELIHOOO FUNCTION AT MLE = –160.29 
SAMPLE FITTED GOODNESS OF FIT TESTS 
MEAN : 3104.4 3104.4 KOLM-SMIRN : 0.20696 
STO OEV : 2944.1 2638.5 SIGNIF : 0.4996 
SKEWNESS: 1.9927 1.6998 CRAMER-V M : 0.08997 
KURTOSIS: 6.8261 7.3342 SIGNIF : > .15 
« BASED ON MIOPOINTS OF FINITE INTERVALS ANOER-DARL : 0.4847 
SIGNIF : "e 
PERCENT ILES SAMPLE FITTEO 
5t 181226 ЖУА KS, AD, AND CV SIGNIF. LEVELS МОТ 
10: 481.26 545.54 EXACT WITH ESTIMATED PARAMETERS. 
291 1259.4 1185.1 
50: 3018 2397.2 NOTE: A SMALL SIGNIFICANCE LEVEL 
253 3344 4267.4 (ЕС. Р<.01) INDICATES LACK OF FIT 
90: 7105.5 5597.3 
95: 12178 8308.3 


Figure 39 Data Table 


ANALYSIS OF WEIBULL OISTRIBUTION FIT 


OATA :  TFUELO8 

SELECT I ON : ALL 

X AXIS LABEL: Dollors 

SAMPLE SIZE : 16 

CENSOR I NG : | NONE 

FREQUENCIES : 1 

EST. METHOD : MAXIMUM LIKELIHOOO 

CONF METHOO : ASYMPTOTIC NORMAL APPROXIMATION 


CONF. INTERVALS COVARIANCE MATRIX OF 
(95 PERCENT) PARAMETER ESTIMATES 
PARAMETER ESTIMATE LOWER UPPER C a 
C (SHAPE) 1.7563 1.0625 2.45 0.12524 8.1824E2 
a (SCALE) 49874 352356 64511 818.24 5.5747E7 
LOG LIKELIHOOO FUNCTION AT MLE - —184.17 
SAMPLE FITTEO GOOONESS OF FIT TESTS 
MEAN : 44424 44409 KOLM-SMIRN : 0.13903 
STD DEV : 26742 26105 SIGNIF : 0.91656 
SKEWNESS : 0.43147 0.9342 CRAMER-V M : 0.032114 
KURTOSIS: 2.3588 1.3818 SIGNIF : > .15 
* BASED ON MIOPOINTS OF FINITE INTERVALS ANDER-OARL : 0.211 
SIGNIF : >415 
PERCENTILES SAMPLE FITTEO 
5: 6565.9 9191.6 KS, AD, ANO CV SIGNIF. LEVELS NOT 
10: 11672 13848 EXACT WITH ESTIMATED PARAMETERS. 
25: 23229 24535 
50: 43152 40480 NOTE: A SMALL SIGNIFICANCE LEVEL 
25: 62744 60068 (EG. Ps.01) INOICATES LACK OF FIT 
90: 89412 80189 
95: 96866 93151 


(21 


DATA 
SELECT | ON 

X AXIS LABEL: 
SAMPLE SIZE : 
CENSOR | NG қ 
FREQUENCIES : 
EST. METHOD : 
CONF METHOO : 
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Figure 43 Data Table 


ANALYSIS OF LOGISTIC DISTRIBUTION FIT 


TSUBSO7 

ALL 

Dollors 

16 

NONE 

1 

MAXIMUM LIKELIHOOD 

ASYMPTOTIC NORMAL APPROX IMAT ION 


CONF. INTERVALS 
(95 PERCENT) 
PARAMETER ESTIMATE LOWER UPPER 
ALPHA 3432.7 2757.3 4108.1 
BETA 795.66 469.56 1121.8 
LOG LIKELIHOOD FUNCTION AT MLE = -140.47 


COVARIANCE MATRIX OF 
PARAMETER ESTIMATES 
ALPHA 
1.176765 e 
0.000Е0 27670 


ВЕТА 


SAMPLE FITTED GOODNESS OF FIT TESTS 
MEAN 3669.8 3432.7 KOLM-SMIRN : 0.23775 
STD OEY : 1894 1443.2 SIGNIF 278: 52626 
SKEWNESS : 1.9387 0 CRAMER-V М: 0.2226 
KURTOSIS: 7.6241 4.2 SIGNIF : > 215 
» BASED ON MIDPOINTS OF FINITE INTERVALS ANDER-DARL : 1.2804 
SIGNIF : DS 
PERCENTILES SAMPLE КҮТТЕ@ 
S: 1039.3 1090 KS, AO, AND CV SIGNIF. LEVELS NOT 
10: 1360; 7 1684.5 EXACT WITH ESTIMATED PARAMETERS. 
Z9: 2113.2 259586 
50: 3985.1 3432.7 NOTE: A SMALL SIGNIFICANCE LEVEL 
75: 3755 4306.8 (EG. Ps.01) INOICATES LACK OF FIT 
90: 5662.2 5181 
95: 9659.5 9779.9 


Figure 47 Data Table 


ANALYSIS OF LOGNORMAL DISTRIBUTION FIT 


OATA TPORTO8 
SELECTION ALL 
X AXIS LABEL: Oollors 
SAMPLE SIZE : 16 
CENSORING : NONE 
FREQUENCIES : 1 
EST. METHOO : MAXIMUM LIKELIHOOO 
CONF METHOO : EXACT 
CONF. INTERVALS COVARIANCE MATRIX OF 
(95 PERCENT) PARAMETER ESTIMATES 
PARAMETER ESTIMATE LOWER UPPER MU SIGMA 
MU 10.253 9.727 10.778 0.057004 O 
SIGMA 0.95502 0.72861 1.5266 0 0.028502 


LOG LIKELIHOOD FUNCTION AT MLE = —186.01 


SAMPLE FITTED GOOONESS OF FIT TESTS 
MEAN 55316 44741 KOLM-SMIRN : 0.28052 
STO OEV : 98294 54603 SIGNIF 0.16115 
SKEWNESS: 3.036 5.4791 CRAMER-V M : 0.26189 
KURTOSIS: 11.176 84.858 SIGNIF : 2:15 
« BASED ON MIDPOINTS OF FINITE INTERVALS ANDER-DARL : 1.3477 
SIGNIF : о 
PERCENT ILES SAMPLE FITTEO 
5: 9.2892E3 5.8922E3 KS, AD, AND CV SIGNIF. LEVELS NOT 
10: 9.6902E3 в. 55 EXACT WITH ESTIMATEO PARAMETERS. 
25: 1.7410E4 1.4894E4 
50: 2.0855Е4 2.8356E4 NOTE: A SMALL SIGNIFICANCE LEVEL 
75: 2.9579E4 5.3985Е4 (EG. Ps.01) INDICATES LACK OF FIT 
90: 1.3051€5 9.6441E4 
95: 4.0299E5 1.3646E5 
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Figure 51 Data Table 


ANALYSIS OF LOGNORMAL OISTRIBUTION FIT 


OATA TEOUIPOS 
SELECT ION ALL 
X AXIS LABEL: Do lors 
SAMPLE SIZE : 16 
CENSOR I NG : NONE 
FREQUENCIES : 1 
EST. METHOD : MAXIMUM LIKELIHOOD 
CONF METHOD : EXACT 
CONF. INTERVALS COVARIANCE MATRIX OF 
(95 PERCENT) PARAMETER ESTIIAATES 
PARAMETER ESTIMATE LOWER | UPPER MU SIGMA 
MU 7.4952 6.7489 8.2416 0.11493 0 
SIGMA 1539551 1.0346 2.1676 0 0.057466 


LOG LIKELIHOOO FUNCTION AT МЕЕ = -147.5 


SAMPLE FITTEO GOOONESS OF FIT TESTS 
MEAN 2682.3 4512.8 KOLM-SMIRN : 0.13769 
STO DEV : 3941.3 10379 SIGNIF 0.92294 
SKEWNESS : 1.0952 19.0656 CRAMER-V M : 0.050035 
KURTOS!S: 2:9511 2178.4 SIGNIF i э 105 
» BASEO ON MIDPOINTS OF FINITE INTERVALS ANDER-DARL : 0.31448 
SIGNIF : ЭТТЭ 
PERCENT IILES SAMPLE FITIEO 
9: 158.08 193.29 KS, AO, ANO CV SIGNIF. LEVELS NOT 
10: 185.64 316.43 EXACT WITH ESTIMATED PARAMETERS. 
291 666.93 721.24 
50: 2541 1799.4 NOTE: A SMALL SIGNIFICANCE LEVEL 
75: 5279.6 4489.3 (EG. P<.01) INDICATES LACK OF FIT 
90: 10486 10233 
95: 12568 16751 
ы - 
Figure 55 Data Table 
ANALYSIS OF LOGNORMAL DISTRIBUTION FIT 
DATA TVREPO6 
SELECT I ON > ALL 
X AXIS LABEL: Oollars 
SAMPLE SIZE 16 
CENSOR ING NONE 
FREQUENCIES = 1 
EST. METHOD : MAXIMUM LIKELIHOOD 
CONF METHOD : EXACT 
CONF. INTERVALS COVARIANCE MATRIX OF 


(95 PERCENT) PARAMETER ESTIMATES 


PARAMETER ESTIMATE LOWER UPPER SIGMA 
MU 11.082 10.495 11.67 0.071287 0 
S IGMA 1.068 0.81479 1.7071 0.035644 
LOG LIKELIHOOD FUNCTION AT MLE = -201.07 
SAMPLE FITTEO GOOONESS OF FIT TESTS 
MEAN 1.1140Е5 1:150ЛЕ5 KOLM-SMIRN : 0.12223 
SID.DEV : 1.1962E5 1.6779Е5 SIGNIF 0.9706 
SKEWNESS:  1.3067E0 7.4826E0 CRAMER-V М : 0.052749 
KURTOSIS: 3.4176E0 1.854262 SIGNIF : prs 
* BASED ON MIDPOINTS OF FINITE INTERVALS ANOER-OARL. : 0.3337 
SIGNIF Pilato 
PERCENTILES SAMPLE FITTEO 
Э: 8.0851E 3 1.1219Е4 KS, AO, AND CV SIGNIF. LEVELS NOT 
10: 2.3280Е4 1.6541Е4 EXACT WITH ESTIMATED PARAMETERS. 
29: 2.6833Е4 3. 1647Е4 
30: 5.4715Е4 6.5019Е4 NOTE: А SMALL SIGNIFICANCE LEVEL 
Z9: 1.4328E5 1::5358E»5 (EG. P<.01) INDICATES LACK OF FIT 
90: 23.1972E5  2.9998E» 
95: 3.9985E5.  3.7881E»5 
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Figure 67 Data Table 


ANALYSIS OF LOGNORMAL DISTRIBUTION FIT 


DATA :  TVREPO6 
SELECTION > АСЕ 

X AXIS LABEL: Dollars 
SAMPLE SIZE : 16 
CENSOR ING : NONE 


FREQUENCIES : 1 
EST. METHOD : MAXIMUM LIKELIHOOD 
CONF METHOD : EXACT 


CONF. INTERVALS COVARIANCE MATRIX OF 


(95 PERCENT) PARAMETER ESTIMATES 
PARAMETER ESTIMATE LOWER UPPER MU SIGMA 
MU 11.082 10.495 11.67 0.071287 0 
SIGMA 1.068 0.81479 1.7071 0 0.035644 
LOG LIKELIHOOD FUNCTION AT MLE — —201.07 
SAMPLE FITTED GOODNESS OF FIT TESTS 
MEAN : 1.1140€5 1. 1501Е5 KOLM-SMIRN : 0.12223 
STO DEV : 1.1962Е5 1.6779Е5 SIGNIF : 0.9706 
SKEWNESS: 1.3067E0 7.4826E0 CRAMER-V M : 0.052749 
KURTOSIS: 3.4176Е0 1.834202 SIGNIF : > .15 
* BASED ON MIDPOINTS OF FINITE INTERVALS ANDER-DARL : 0.3337 
SIGNIF : > .15 
PERCENTILES SAMPLE FITTED 
5: 8.0831E3 1.121964 KS, AD, AND CV SIGNIF. LEVELS NOT 
10: 2.528064 1.654164 EXACT WITH ESTIMATED PARAMETERS. 
25: 2.683364 3.164764 
50: 5.4715E4 6.501964 NOTE: A SMALL SIGNIFICANCE LEVEL 
75: 1.4328Е5 1.3358Е5 (ЕС. Р<.01) INDICATES LACK OF FIT 
90: 3.1072E5 2.5558Е5 
95: 3.9986E5 3.7681Е5 


Figure 68 Data Table 


ANALYSIS OF LOGNORMAL DISTRIBUTION FIT 


DATA : FULLVAL2[;7} 
SELECTION : ALL 

X AXIS LABEL: Dollors 
SAMPLE SIZE : 29 
CENSORING : NONE 


FREQUENCIES ; 1 
EST. METHOO : MAXIMUM LIKELIHOOD 
CONF METHOO : EXACT 


CONF. INTERVALS COVARIANCE MATRIX OF 
(95 PERCENT) PARAMETER ESTIMATES 


PARAMETER ESTIMATE LOWER UPPER MJ 5 ICMA 
MU 11.236 10.868 11.604 0.03115 0 
SIGMA 0.95045 0.7676 1.30B3 0 0.015575 
LOG LIKELIHOOD FUNCTION AT NLE = -365.51 
SAMPLE FITTED GOOONESS OF FIT TESTS 
MEAN — : 1.1538Е5 1.1904Е5 CHI-SQUARE : 1.8769 
STD DEV : 1.0854Е5 1.4423E5 DEG FREED: 1 
SKEWNESS: 1.3233Е0 5.4130E0 SIGNIF  : 0.17069 
KURTOSIS: 3.6926Е0 8.242261 KOLM-SMIRN : 0.095754 
* BASED ON MIDPOINTS OF FINITE INTERVALS SIGNIF  : 0.95304 
CRAMER-V M : 0.042478 
PERCENTILES SAMPLE FITTED SIGNIF : > .15 
5: 2.528084 1.586564 ANDER-DARL : 0.28923 
10: 2.354584 2.241264 SIGNIF : > .15 
25: 3.755664 3.992664 
50: 6.615064 7.577864 KS, AD, ANO CV SIGNIF. LEVELS NOT 
75: 1.5418ES 1.4382Е5 EXACT WITH ESTIMATED PARAMETERS. 
90: 3.1072E5 2.5622Е5 
95: 3.6854E5 — 3.6194E5 NOTE: A SMALL SIGNIFICANCE LEVEL 
(EG. P£.01) INDICATES LACK OF FIT 
CHI-SQUARE GOOONESS OF FIT TABLE 
LOWER UPPER 085 EXP oŒ . ((0-E)*2)-E 
-INF. 6.5297E4 14 12.695 1.3047 0.13408 
6.5297Е4 1.305985 5 8.085 -3.085 1.1771 
1.3059Е5 1.958965 5 3.6135 1.3865 0.53198 
1.958985 +ІМЕ. 5 4.6061 0.39386 0.033677 
TOTAL 29 29 1.8769 
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' Edit Time Assumption Dialog Box Module 


" LCDR Terry Redman, USN 
' Naval Postgraduate School 
2994 


ХХХ ЕЕЕ ЕКЕ eode oe ode oe oe oe eoe ode eoe obe eoe oe de eoe ode oe oe oe eoe eoe oe de Ж Ж Ж Жз Жусу 


puttoni9 Click Macro 


Worksheet Edit Button Actions 


EupEEuttonl9 Click() 


DialogSheets("Dialogl").Show 
DralogSheets("Dialogl").OptionButtons(" "Option Button 20") - xlOn 
Prabogsnecets (“Dialogl") .EditBoxes ("Edit Box 23") .Text = "1" 
DialogSheets ("Dialog1") .EditBoxes ("Edit Box 14").Text = " " 

Maa togoheets ('Dialogl”).EditBoxes ("Edit Box 19").Text - 


= 
-. 
= 
- 


Ena Sub 


к. ХЕ ККЕ e obe e eoe e oe ode oe Жу eoe oe oe oe fe oe oe ode oe oe ode ode oe oe oe oe ode e de oe oe e obe eoe de oe ok Kh +) Keke Ж Ж Ж Ж Ж Ж Ж б Ж АЕ 


Poeton Click Macro 


OFS BUttOon actions 


1 
1 
1 
1 
! 
1 


EnpSPuUtton2 Click() 


If DialogSheets("Dialogl").OptionButtons("Option Button 20").Value - xlOn Then W 
orksheets("Sheetl").Cells(36, 3).Value - DialogSheets("Dialogl").EditBoxes("Edit Box 
234 Text 

If DialogSheets("Dialogl").OptionButtons("Option Button 11").Value = xlOn Then W 
orksheets("Sheetl").Cells(36, 3).Value - ((DialogSheets("Dialogl").EditBoxes("Edit B 
Text) * (0.032967033)) 

If DialogSheets("Dialogl").OptionButtons ("Option Button 15").Value - xlOn Then И 
orksheets("Sheetl1").Cells(36, 3).Value - ((DialogSheets("Dialogl").EditBoxes("Edit B 
РЕ ТехЕ) * (0.2307692308)) 


Ema Sub 


ХЕ ХЕ ЕЕЕ ЕЕЕ ЖЕКЕ ЕЕЕ ЕЕЕ eode kkk KKK kK Ж} 


i Button Click Macro 


1 
' Cancel button msg 
EubSButton$ Click() 
MsgBox "Are you sure you Want to Cancel?", vbYesNo 


End Sub 


Y eoe e e je e e e de de e e de ke e e e e e e e e e de e e oe e oe oe e ke e e oe ck oe de oe oe de ode oe e de e ke ose ke oe de oe e e de de de oe dee eoe oe oe ode obedece ck de eee Б АЕ 


" EditBox14 Change Macro 


ШГ. < Edit Box Actions 
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Sub EditBoxl4 Change () 


DialogSheets ("Dialogl").OptionButtons ("Option Button 11") = х10п 
DialogSheets ("Dialogl") .OptionButtons ("Option Button 20") = xloOff 
DialogSheets ("Dialogl") -OptionButtons ( Opticn Button у = ое 


MsgBox "Are you sure you want to make this change?", vbYesNo 


End Sub 


Te ke oe ce oe oe oe ck oe e ce oce eoe che ode oe ode oe К А А А А А А К А А eoe ode ose eode oce ode oe oko ok oe oe oce ok ook ook oe ook ode oe ode ЗАС АС АС С С АС АС АС oe oce ode с С АС ko kx kx AX 
1 


' БкаївВох19 Change Масго 


1 


' Weeks Edit Box Actions 


! 


Sub EdrtBoxlo рана 


DialogSheets("Dialogl").OptionButtons("Option Button 15") - xlOn 
DialogSheets ("Dialog1").OptionButtons ("Option Button 11") = xlOff 
DialogSheets ("Dialog1™) .OptionButtons ("Opt oDFBUECODTZ2 OTME em 


MsgBox "Are you sure you want to make this change?", vbYesNo 


End Sub 


Y oe ke ode oce ode oe oce ode oe oe oec oce oce oe oe ode ok oe ode oe ok eoe dee ke ck ode oec oe oe oe oe ode oe ode ode ck eoe oe ode ek ck oe oe ck oe oe oko +» +» +)” +» +» К А К А К Ж А А Ж А € 


' EditBox23 Change Macro 


1 


t Months Edit Box Actions 


1 


Sub EditBox23 Change() 


DialogSheets("Dialogl"™) .OptieonButtons @ Option Buttons О Е Сом 
DialogSheets ("Dialog1").OptionBüttons ("Option Button I) = Бои 
DialogSheets ("Dialog1") .OpticonButtoóons (“Option Butten 15.) — ier 


MsgBox "Are you sure you want to make this change?", vbYesNo 


End Sub 


ТХХХХХхХхХХХххХХ eoe oe oe ck ck oe ke ke +» +» +) +» oce oko oe К А К А А А Ж К А К К К ХХХ Хе жж 


1 


' OptionButtonll CIickiMacre 


1 


' Days Option Button Action 


1 


Sub OptionButton]l] Clem 


DialogSheets ("Dialogl") -OptionEuttons ("Cptreonmeurron. 20) xLOLG 
DialogSheets ("Dialogi") OptionBüttons| ODESODEBUEESODEMISUDILCEPTE 


End Sub 


к КЕК ЕКЕ + АС +» С С с А А +) +) +) +» б е б А К А А К С А б А К С В ЯС С АЕ К + + + + Ж ЖК 


' 


' OptionButtonl> ClickeMacre 


! 


' Weeks Option Button Action 


sub Opti onButton) > Се 
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DialogSheets("Dialogl").OptionButtons("Option Button 20") XJO0ff 
REnsTogsShestsi"Dialogl" ).OptrionButtons("Optazon Button 1l") 9 xlofft 


End Sub 


e e se e e Je de e e e de Je e ode oe oe oe oe oe oe oe eode ode oe de oe ode oe ode oe ose ode oe ode eoe oe oe oe eoe oe eode eode oe eode ode АС eode eoe oe oe с Эс eoe eoe eoe Ус Ус eo eo ЭЖ Ж ЗУ ЖК 


mogetonButtonz0) Click Macro 
! 
 Мешенс Option Button Action 


Sub OptionButton20 Click() 


pudosgsheets("Dialogl").OptionButtons( Option Button 11") OPI 
DialogSheets ("Dialogl") .OptionButtons ("Option Button 15") = xloff 


End Sub 
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' Run Dialog Boz Module 

' LCDR Terry Redman, USN 

' Naval Postgraduate School 
' 1994 

"aButton299c]gc асса 

t 

' Worksheet Button Actions 
Subi Büttonz2? CrrEkI) 


DialogSheets ("Dialog4") .Show 
DialogSheets("Dialog4").EditBoxes("Edit Box 6").Text — "2000" 


End Sub 


Buttons likek Macro 


OK Button Actions (Reserved) 


Sub Button? Crrek() 


End Sub 


ВЕС Macro 


Cancel button msg 
Sub ButtonoSwc rci 
MsgBox "Are you sure you want to Cancel?", vbYesNo 


End Sub 
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' Edit Budgeted Overhead Dialog Box Module 


' LCDR Terry Redman, USN 
' Naval Postgraduate School 
1474 


«Je Je de e KK a Е Уе Ж Ус К % 


! 


: Breton Click Macro 


t 


' Worksheet Button Actions 


Sepeeucton2/ Click() 


DialogSheets ("Dialog2").Show 
EEogSheets("Dialog2" ).OptionButtons("Option Button 5").Value - xlon 


DialogSheets("Dialog2").EditBoxes("Edit Box 8").Text = "0" 
wPalogsheets("Dialog2").EditBoxes("Edit Box 12").Text = "0" 
End Sub 


ЕКЕ Уб e eoe de eode de de ode oe oe oe oe ode eoe oe oe e eoe ode de oe eoe eoe de Ke дж 


! 


І Poe on Click Macro 


' OK Button Actions (Reserved) 


EubOButton2 Click() 


ibrepialogSheets("Dialog2").OptionButtons("Option Button 5").Value = xlOn Then Wo 
Bksheets("Sheetl").Cells(32, 3).Value = (((DialogSheets("Dialog2").EditBoxes("Edit B 
ox 8").Text / 100) + 1) * Worksheets("Sheetli").Cells(32, 3).Value) 

If DialogSheets("Dialog2").OptionButtons ("Option Button 9") .Уа1ае = xlOn Then Wo 
rksheets("Sheetl").Cells(32, 3).Value = ((1 - (DialogSheets("Dialog2") .EditBoxes ("Ed 
it Box 12").Text / 100)) * Worksheets ("Sheet1").Cells(32, 3) .Value) 


End Sub 


V Je ke e Je dee Ееее Ееее ж 


! 


{ Botton Click Macro 


' Cancel button msg 
1 


suo Buttons Click() 
Mseqbox “Are you sure you want to Cancel?", vbYesNo 


End Sub 


Y Je Je Je e de Je e de dee oe e de ke oe de e de e ode eode ke e dede e de oe de ЕЕЕ Ее еәеәж.ж 


BeceeronButtons Click Macro 


Sub OptionButton5 Click(í) 
prualogSheets('" Dialog2").OptionButtons("Option Button 9").Value - xlOff 


End Sub 
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Vox oe e de de e e he ode de oe ce ode he oe eoe e dee ode oe oce eode he ode oe oe eode eode oe kk k k de ode oe ode oe hee ke oe eode ode oe ode eoe eoe + + k k k ke k k k ke ke ke ke ke ke ke +++) Ж 


' OptzonButcond Glick Macro 


SuboOptionbuccon mel pcr 


DialogSheets ("Dialog2").OptionButtons("Option Button 5").Value - xlOff 


End Sub 
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' Edit Number Of Ships Dialog Box Module 


' LCDR Terry Redman, USN 
' Naval Postgraduate School 
P 1994 


foe oe Jede oe ode oe ode oe oe ode eode ode ke oe ode oe oe oe e ode ode eode de ode ode ode ode oe e ek oe ode oe koe oe oe oe oe oe eode eode oe eode Ж Ж * * Ж К Ж Ж + Ж Ж Ж * Ж de Ж Ж Ж} Ж Ж * Ж АЕ 


: Buticemz © Click Macro 


' Worksheet Button Actions 


' 


Sub Butcon2” Click (} 


DialogSheets ("Dialog3") .Show 
DialogSheets ("Dialog3") .EditBoxes ("Edit Box 6").Text = "1" 


End Sub 


Y eoe ode de oe de ode oe oe oe ode e ЕЕ ЕХ АА << еее ж 


Eaelog5 Button2 Click Macro 


' OK Button Actions 


! 


ООШ акбоз Buttonz Click() 


Worksheets("Sheeti1").Cells(39, 3).Value = DialogSheets ("Dialog3") .EditBoxes ( 
ШЕСІШЕ Бок 6”).Текхі 


End Sub 


Y exe e de de de de de ke de de e de e oe koe e e oe oe de e oe ЕЕЕ Хе ЕЕ ЕЕЕ Хаж 


| Bncbonsoclick Macro 


! 
E cancel button msg 
BubcButton3 Clicktí) 


MsgBox "Are you sure you want to Cancel?", vbYesNo 


End Sub 
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' Crystal Ball Macro Module 


' LCDR Terry Redman, USN 
' Naval Postgraduate School 
' 1994 


Toe xe ode de e eoe ode oe oe ode ode eoe ode ode eode ode oe oe eoe ode eoe ode Ж eoe ode oe ode ode ode oe Ж 5 Ж Ж + Ж кх Ж ЖЖ Ж з + + Ж + Ж eode oe oce eode oe + Ж + Ж ХХ Ж Ж + 


' AccessAssumption Macro 
' Macro recorded 4/10/94 by Terry Redman 


1 


Sub AccessAssumption() 
Application.Run Macro:=Range("CB.DefineAssum") 


End Sub 


1 Ж ke ook к кке е 


' CreateReport Macro 
' Macro recorded 4/10/94 by Terry Redman 


Sub CreateReport() 
Application.Run Macro:-Range("CB.CreateRpt") 


End Sub 


ХХ + Ж К К К Ж К Ж К ж Ж ж К ж ж К Ж Ж + Ж Ж Ж + Ж Ж Ж Ж Ж ++ + Ж} Ж Ж + + К oe Ж + + Ж + + + Ж + + + ЖЖ Ж 


' RunReplications Macro 
' Macro recorded 4/10/94 by Terry Redman 


Sub RunReplications() 
Application.Run Macro:-Range("CB.Run") 


End Sub 


ЕЕЕ е е ЕКЕ aK 


" ResetRun Macro 
' Macro recorded 4/10/94 by Terry Redman 


Sub ResetRun() 
Application Run Macro; —Range( Gb Reser | 


End Sub 
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